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ABSTRACT 


This work deals with the evolution of a FORTRAN 
Simulation written by Naval Research Laboratory which is 
used to evaluate the effectiveness of deploying ship- 
launched infrared decoys to counter the anti-ship infrared 
seeking missile threat. Although the model (referred to as 
the Stochastic Infrared Engagement Model —- SIREM) possesses 
extensive analytical capability and flexibility, refinements 
are desired to more accurately emulate atmospheric effects 
on the acquisition process. Methods are derived herein to 
calculate atmospheric transmittance as a function of range 
using an accurate, LOWTRAN-based empirical formula. 

Basic seeker discrimination techniques are addressed 
which may be incorporated into SIREM or other simulations 
for future missile-versus-decoy evaluations; and some ideas 
are presented which may prove valuable in decoy enhancement 


to subjugate the discriminating seeker. 
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I. INTRODUCTION 

With the advent of the missile era has come a revolution 
in naval warfare methods. The problem of contending with 
"smart" weapons has given a whole new meaning to both 
offense and defense alike. These missiles may employ a 
variety of acquisition hardware. Some will be passive, some 
active, and some both. 

Farticular problems arise when missiles use passive 
acquisition methods. Passive sensors radiate nothing and 
are therefore invisible to conventional electronic support 
measures (ESM). Furthermore, their small size allows them 
to reflect only meager amounts of radar energy, which 
further reduces reaction and countermeasure employment time 
of the ship. The countermeasures may be in the form of 
onboard or offboard weapons systems or deception devices. 
The cost versus effectiveness of onboard close-in weapon 
systems has been a driving factor in the Navy' s desire to 
analyze the employment of offboard expendable seduction and 
distraction devices as a relatively inexpensive alternative. 

Here we will restrict our studies to the passive 
infrared seeking missile problem, although concurrent 
analysis is being conducted to contend with the active RF 


missile seekers. 


To reduce radically the cost of development of offboard 
device characteristics and deployment tactics, a computer 
simulation has been constructed to assist in the analysis 
process. This modelling program is appropriately named the 
"Stochastic Infrared Engagement Model - SIREM". It 
originated at Naval Kesearch Laboratory in Washington, D.C., 
and was transported to the Naval Postgraduate School in 
Monterey, CA to assist in the analysis of the use of 
infrared decoys as a counter to the passive seeker problem. 
It is designed to imitate the typical cruise-missile versus 
ship engagement scenario and to establish measures of 
effectiveness (MOES) for both 1) decoy characteristics, and 
2) consue tactics. Therein lies the basis for 
continuing verification and upgrades to the model so that it 
remains capable of accurate evaluation of the hardware at 
hand. 

Most real-world systems are very complex. In order to 
model any system accurately, all factors which affect it 
must be represented in as accurate a fashion as possible. 
Weaknesses in representation should be identified and 
documented so that overconfidences in the output do not 
occur. Once documented, the weak areas of the model may be 
studied and reworked to increase the accuracy of the model. 

This paper deals with the evolution of SIREM. Εν 
Simulating each component of the engagement scenario in a 


realistic manner, the optimum decoy deployment schemes can 


be obtained, combined with their RF counterparts, and 
algorithmically automated for deployment via computer in 
order to reduce chances of operator control errors. It can 
easily be seen here that it is not financially or 
logistically feasible to evaluate all scenarios in field 
tests using real ships and missiles, since the cost would be 
prohibitive. 

Extensive conversion and review of the SIKEM model was 
conducted to assist in its enhancement. Two specific areas 
Of weakness were identified which require some attention. 
The first area of concern was a lack of atmospheric 
attenuation on the emitted IR signatures of the ship and/or 
decoys. The second regards the study of seeker head 
discrimination techniques and how to implement them based on 
data available. The emphasis in this paper revolves around 
finding a simple technique for calculating atmospheric 
transmittance throughout the path of an incoming infrared 
seeking missile which will accurately affect its probability 
of acquisition in the simulation environment. Two methods 
currently used in obtaining transmittance in a stochastic 
model such as SIREM involve either look-up tables which are 
either incomplete or inflexible, or incorporation of the 
LOWTRAN program as a subroutine. The second method is 
preferred since it is very flexible to observed atmospheric 
conditions; however, it requires an extension of 9000 lines 


of FORTRAN code along with the added input-output coding, 


which serves only to calculate one value of atmospheric 
transmittance for each incremental range throughout the 
missile flight path. It was felt by the author that this 
was excessively wasteful of valuable computer time and 
memory requirements considering the model's current size of 
approximately 15,000 lines of code. 

The intent was to devise a simpler scheme for finding as 
accurately as possible what the atmospheric transmittance 
would be for a given set of atmospheric conditions that 
would be easily incorporated into the incremental flight 
path of a missile and more accurately affect the seeker head 
acquisition probability. 

To more fully understand the principle of atmospheric 
transmittance and its effects, one chapter is dedicated to 
explaining the basic principles of atmospheric 
transmittance. Once a good understanding of the principles 
of atmospheric transmittance has been gained by the reader, 
an introduction to the LOWTRAN program is performed, 
presenting its capabilities for prediction of transmittance. 
This will be followed by a simplistic approach for modeling 
atmospheric transmittance with very short computer coding 
schemes. 

The closing chapters will be dedicated to more detailed 
explanation of the SIREM modelling program, followed by an 
introduction to possible discrimination techniques Since 


seeker discrimination adds to the problems surrouding the 


use of offboard countermeasures, the subject of enhancing 
current decoy designs by taking advantage of certain natural 
phenomenon is addressed in the concluding chapter. 

Future enhancements to SIREM to incorporate both 
transmittance and discrimination techniques are the object 
of this work and concurrent work by Naval Research 
personnel. These enhancements are intended to more 
accurately reflect the acquisition probability of the 
seekers being modelled so that refinements to the design and 


deployment of offboard countermeasures are possible. 
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II. ATMOSPHERIC ENVIRONMENT 

This chapter is dedicated to reviewing the principles of 
atmospheric propagation of infrared radiation, the causes of 
its variability, and tools available to predict how much is 
transmitted at various ranges and atmospheric conditions. 

The atmosphere is comprised of numerous gasses - 
nitrogen, oxygen, water vapor, carbon dioxide, methane, 
carbon monoxide, and ozone to name the major components. 
Additionally there are numerous suspended particles 
distributed in various densities throughout the earth s 
atmosphere. The major constituents of the atmosphere are 
nitrogen, comprising 784, and oxygen, which occupies 20% of 
the total volume. Water vapor makes up 17, while carbon 
dioxide constitutes only 0.044 of the atmosphere. The 
transmission of infrared radiation through the atmosphere 
depends on the meteorological conditions, and hence varies 
with weather conditions and altitude. Only an approximate 
result can be gained through the theoretical study of the 
Principles of atmospheric attenuation since no definitive 
conclusions can be made. The base of knowledge regarding 
atmospheric attenuation and the associated formulae comes 


principally from empirical means or approximations. 


TL 


The primary factor in absorption in the infrared 
spectrum is water vapor. The wide variations in 
concentration and distribution of water vapor in the 
atmosphere are the primary reasons for the variance in 
absorption levels from one meteorological profile to the 
next. Although carbon dioxide comprises only 4/100 percent 
volume of the atmosphere, it is second in importance as an 
attenuator of infrared radiation, as well as being more 
uniformly distributed than water vapor. 

There are three major factors which contribute to the 


attenuation of infrared radiation: Extinction, or the loss 


on! etre À 





of energy resulting from the interaction of the beam with 
the various absorbing or scattering constituents contained 
in the atmosphere; Refraction, or the bending of light rays 
due to the refractive index gradients of the atmosphere, and 
Scintillation, or the distortion of the optical beam due to 
small scale turbulence. Refraction will not be discussed 
heavily here since it will be apparent later that it may 
serve only to extend the optical horizon of infrared seeking 
missiles, but have little affect on the simulation discussed 
herein. Extinction represents the combined effects of the 
scattering of radiation into or out of the beam due to 
aerosols and air molecules present in the beam path, and 


absorption caused by the direct transfer of energy from 


radiation generated by the source to the vibrational or 


o 
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rotational kinetic energy of the molecules and aerosols. 


This is illustrated in the following diagram. 
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Figure 2-1. Causes of extinction 


The attenuation or extinction of electromagnetic 


radiation is described by the Lambert-Beer-Bouguer Law 
f= [ores (EG. 2) 


where I is the attenuated radiation, I, 15 the source 
radiation, f is the attenuation coefficient, and KH is the 
path length. The extinction is caused by absorption and 
Scattering by molecules and aerosols. They contribute 


linearly to the the total extinction as follows: 


bh = bua + Pns + Faa + Bas (Eq. 2.2) 
where Z4 molecular absorption by water vapor, ozone, etc. 
βκεξ molecular (Rayleigh) scattering 
Da aerosol absorption (dry particles neglected) 


Pas aerosol (Mie) scattering 


The importance of each of these effects depends on the 
wavelength of the electrooptic energy. The wavelength 


ν 
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dependencies exist because of the relative sizes of the 


responsible constituents (molecule or aerosol) and the 


energy states of the molecules (molecular absorption bands). 


The importance of the principal component of attenuation for 


the primary wave bands is given here. 


TABLE 2.1. PRIMARY ATTENUATORS IN IR SPECTRUM 


Wavelength Region Attenuation Coefficients 

_( In micrometers ) In order of importance 
Visible .4 — .7 micron Pra, Ons 

Near IR .7 - 2 micron Bas (Ema for some wavelengths) 
Mid IR 5 = micron 17৭50171250) 5 Bas 

Far IR 8 - 14 micron fun (FRO), Bas 


The aerosol relationship used to describe extinction is 


Principally governed by the size (aerosol radius) of the 


Particle and the wavelength of the incident energy as: 


1) Rayleigh Scattering (where r<<} — molecule effects). 
2) Mie Scattering (where r £$ > approximately). 


3) Geometric Scattering (where r»^^À — large particles). 


An illustration of the scattering description is as shown 


here. 
[($)2Scattered light 






Scattering hngle 


I = Incident Energy | 9=0° 
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Figure 2-2. Scattering Geometry 


14 


The attenuation scattering cross-section for one particle of. 
size r and wavelength à is the integral of the scattering in 


all directions: 


4 


5 = 709১ dQ (Ed 10 


Solid Angle 


) 


ہے 


Mie scattering is generally concentrated in the forward 
direction as follows. 


1 (9) 
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Figure 2-3. Mie Scattering 


Note that a minimum occurs near the +100° points and that a 
secondary peak occurs in the backward direction (back- 
scatter). 


Rayleigh scattering is symetrical as shown. 


Figure 2-4. Rayleigh Scattering 


The scattering efficiency Q of a particle is given by the 


cross-section S (above) divided by the area of the particle. 


ar? (Eq. 2.4) 


For Rayleigh scattering Q is on the order αἵ λ΄. and based 


on the wavelength span present from the sun' s radiation, is 
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the principal factor in the sky's blue appearance. AQ is 
higher at the longer wavelengths. For geometric scatterinc 
Qs 2.0. Only particles with radii between about 0.1 and 
1.0 micrometers are effective in scatter’*:iig light in the 
visible range. These particles are too small to rain out 
due to their low colle-cion efficiency. This is why rain 
alone does not rezZuce haze and increase visibility. A 
char.,e in the vertical temperature profile (such as the 
breaking up of an inversion layer by a frontal passage) does 
increase horizontal visibility through upward diffusion of 
the aerosols. 

Because the scattering efficiency Q is a function of 
aerosol size and wavelength, scattering itself is a function 
Of wavelength. 


Pas ۂ‎ AS (Eq. 2:9» 


Here « is the Angstrom coefficient, and varies generally 
from 0.9 to 2-0: This value may be as high as 4.0 for 
Rayleigh scattering. The scattering coefficients for 


various aerosol types are shown below. 
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Figure 2-5 Scattering Coefficient 
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With this diagram we see that the wavelength dependence of 
scattering explains why the sky is blue, since the shorter 
wavelengths are scattered so much more. Likewise heavy haze 
can cause the red sunset effect since the longer wavelengths 
are not scattered out of the path. Note also how the cloud, 
fog, and haze coefficients carry over into the infrared 
region. We can expect these to be major factors in 
transmittance reduction there. 

Scattering is also indirectly a function of relative 
humidity (RH) because as HH increases, the water molecules 
collect on a dry particle until enough of them are on the 
particle to cause it to go into a solution and form a 
droplet. This 1s called deliquescence. Since a solution 
droplet can grow in size much faster than a particle, the 
scattering increases rapidly with RH due to an aerosol 
rapidly changing its size distribution through growth of 
droplets. Typical scattering coefficient variations are 
shown below. These are presented as a ratio of the 
scattering coefficient for aerosols to the scattering 
coefficients for aerosols at 20% relative humidity as a 


function of relative humidity. 


4 
15৪, ১ (5208) 2 
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Figure 2-6 Scattering Coefficient Variation 
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Models of the aerosol distribution in the marine 
environment can be very complex because of the amount and 
size of the particles present in the marine environment. 
These are functions of particle generation at the sea 
surface and the vertical transport of the particles once 
they are suspended in the atmosphere. 

Wind speed dependence on aerosol distribution in the 
maritime model is similar to that of continental models 
until the wind speed nears 14 knots. At 14 knots the 
aerosol concentration has increased by a factor of about 1.5 
due to aerosol generation by the surface whitecaps which 
form at this wind speed. As the wind speed increases, the 
concentration may rise to a value of 10 to 100 times the 
calm atmospheric value near the surface, but decrease in 
concentration with altitude up to approximately 3 
kilometers. 

It can now be seen that the atmospheric conditions have 
a dramatic effect upon the extinction and absorption at 
various wavelengths. In this paper the interest lies 
primarily in the infrared region since this is where our 
model (SIREM) is designed to simulate target acquisitions. 
Hence it is good practice to become familiar with those 
factors which affect the systems being modeled. 
Incorporation of "real-world" occurrences is a vital part of 
Systems modeling. Those models which use random (non- 


deterministic) methods to emulate a patterned parameter can 
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only be termed less than valid. They would serve useful for. 
nothing more than gaming exercises since the answers could 
not be relied upon. 

Since the SIREM model is designed to obtain measures of 
effectiveness of IR decoy deployment versus an IR seeking 
missile, those measures may only be obtained through 
realistic incorporation of the factors that affect the 
engagement. The primary intent of the following chapters is 
to define the scenario, present the principles of detection, 
and show how the LOWTRAN program can be used to develop a 
very compact and accurate method of incorporating 


atmospheric absorption into the acquisition process. 


L9 


III. ENGAGEMENT SCENARI 


As with any detection method, whether RF, visible, or 
infrared, it is assumed that the target is either emitting 
or reflecting energy at some power density in the frequency 


band of the detector as in the scenario shown below. 
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Figure 3-1. IR Missile Attack Scenario 





From extended ranges this emission or reflection can be 
assumed to appear as a point source, however at closer 
ranges the physical size of the target must be expanded into 
a distributed source emitter, especially considering the 
typical optical magnification performed by the seeker. 
Since this document deals specifically with infrared 
emission and transmission, the case of target reflection 
will henceforth be disregarded. 

The unit of measure given to the power emitted from the 
target is the radiant intensity (J), commonly stated in 
watts per steradian (solid angle). The target, assumed to 


be a point source, emits a source level radiation in 
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specific spectral bands. The radiant flux density of any 
wavelength across the continuum has a natural reduction in 
value based on the inverse square law with range from the 
ship (target) to the missile (receiving platform). 

The detector in the seeker head is composed of a 
material which electrically (via photon energy) sensitive to 
incoming radiation across some defined frequency bandwidth 
inherent to the composition. This bandwidth is preferrably 
fairly narrow. If it is allowed to become excessively wide, 
the signal-to-noise ratio is reduced an unacceptable level 
of false alarm detections as discussed by Hudson CRef. 2:ch. 
pz]. For this reason most detector systems have external 
spectral filters to narrow the bandwidth and eliminate these 
misdetections. 

Based on the composition of the detector material each 
detector will have an optimum electrically active wavelength 
to which it is sensitive. As the incoming energy shifts in 
frequency to either side of the optimum, the electrical 
activity in the detector is reduced, creating some 
distribution which reflects the electrical activity, or 
molecular transition density as a function of wavelength. 
This relationship should be known in order to reflect any 
spikes in the sensitivity curve where certain energy bands 


may be particularly absorptive of incident energy. The 
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profile shown in Figure 3-2 1s one which could be used in 
the typical simulation environment based on the above 


principles. 


Detector 

Radiation 
Detectivity 
(Relative) 


2.0 5-0 4.0 .. 


Figure 3-2.Detectivity over Bandwidth 


The actual electrical bandwidth of a detector has three 
Primary dependencies: first, the material from which it is 
constructed; second, the temperature at which it is 
operated; third, the physical area of the detector surface. 
The detectivity (D*) is the convenient measure categorizing 
a detector s operation within these parameters. It 
normalizes the "sensitivity" of various detectors to an 
electrical bandwidth of 1 Hz, and an area of 1 square 
centimeter, and is a convenient means for comparing the 
equivalent electrical response between detectors in circuits 
of different bandwidth. Shown in Figure 3-3 18s a typical 
Lead Sulfide detector's D* characteristics taken from Hudson 


(Ret. 2:pp. 3651. 


22 


Lead Sulfide (PbS) 
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Figure 3-3. Detectivity versus Wavelength 


The value of D" is dependent upon the temperature, 
wavelength and bias current applied. The low temperatures 
are required, as well as bias current, to maximize the 
System's signal-to-noise ratio. 

The object of our simulation is to characterize the 
search and acquisition process as closely to actual as 
possible. To do this we must model the detection process 
(seeker), atmospheric effects, flight characteristics of the 
missile, and target (ship or decoy) signatures as closely as 
possible. Taking all factors into account, the maximum 
range of the detector must be calculated. 

The spectral irradiance from a target (assumed distant 
in the field of view (FOV) is 


H, =J Tall in Watts/cm? (Eq. 3.1) 
R2 
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where 7,4) 15 the transmittance of the path, and the 
subject of following chapters. Ris the distance from the 
seeker to the target. 


The spectral radiant power into the detector is 


F,  ΠιβοΤοί(λ) in Watts (Eq. 3. 


where Ao is the area of the entrance aperature of the 
optics, and 7,.(4) is the transmittance of the whole sensor 
system. Το(λ) is the product of the transfer functions for 
the detector and all external filter functions such as 
protective windows, lenses, reticles, and maskings. 


The signal voltage from the detector is 


Vs — 7১৪ (১) (Eq. 3.3) 


where &()) is the spectral responsivity of the detector. 
Thus far the equations represent an infinitesimal 
spectral interval centered about some wavelength }. Since 
there is a finite bandwidth to an actual system, this should 
be taken into account as shown in Hudson LRef. 2sch. 1351. 
Remember, however, that we are assuming "worst case" in the 
Simulation process. So, by selecting the optimum detectivity 
wavelength of the particular detector, we have effectively 
maximized the range of detection for our seeker. This 
assumption may only be made if the seeker system being 


modelled has a relatively narrow spectral response due to 
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external filtering for a better signal-to-noise ratio. This. 
is virtually always the case from design inception since the 
false alarm problem is critical here. 

Further development of the equations above is performed 


by Hudson and leads to the Optical Trade-Off Equation: 


R = [JTa J'2 سس‎ mae 1/2 [107 32 1. Ν 1/2 
2 (WAT) ‘2 (Vs/ Vn) 
آ+9.‎ N 2 Z Neo ১, 4 / (Eq. 3.4) 
1 — Target irradiance and atmospheric transmittance 
2 — Optics parameters 
5 — Detector characteristics 
4 — System parameters and signal processing 
where R = detection range. 
D. = diameter of entrance of aperture optics. 
NA = numerical aperture [Ref. 2:ch. 51. 
w = instantaneous field of view of sensor. 


Af = frequency bandwidth in Hz. 
Vs/ Vn = signal to noise ratio (Ref. 2:pp. 419]. 

By breaking the equation into sections as shown above we 
see in the first term that reduction of target signature 
(ship) and enhancement of the decoy signature inband would 
shift the acquisition probability more toward the decoy. In 
the engagement process nothing can be done about the 
atmospheric transmittance, although it does have dramatic 
effects on the acquisition range. Knowing the parameters of 
the second, third and fourth term facilitates modelling 


various missile/seeker systems. 
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At ranges beyond the sensitivity of the seeker head, 
i.e. where incident energy is insufficient, detection 
probability would be zero. In addition to the inverse 
square law power loss, the atmospheric losses previously 
discussed must be taken into account in order to accurately 
determine detection probability. Since the detectivity of 
the seeker is independent of the atmospheric effects on the 
emitted signature, the resultant probability of detection 
would be the product of probability of detection assuming no 
atmospheric effects, and the atmospheric transmittance which 
varies from © to 1. 

Few models if any actually emulate the signal 
integrating principle upon which most detectors depend for 
threshold detection. The scanning technique, whether serial 
or parallel, performs the temporal integral of the incoming 
energy across the "window function" which defines the 
acceptance bandwidth of the seeker. The seeker system 5 
probability of detection curve-versus-wavelength is 
proportional to the detector s detectivity-versus-wavelength 
cur ve. 

The present version of the SIREM model emulates a 
scanner type seeker with limited field of view typical of 
the ASMs previously discussed. This model is presently 
being upgraded to account for atmospheric losses over range. 
It will do so by assuming the "worst case" posture. Hy 


this, it is meant that at the optimum detector wavelength, 
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the model presently calculates probability of detection for 
incremental missile movements. This value is calculated as 
a "worst case" value for the ship based on the detectivity 
of the seeker and the maneuvering capability of the missile. 
By multiplying the probability of detection (assuming no 
losses) at each range increment by the corresponding 
transmittance value at that range, the target acquisition 
probability is more accurately represented. This is the 
case whether the acquired target 1s a decoy or the ship. 

As an example assume that the transmittance has been 
calculated for a given atmospheric profile for the 
engagement scenario. The transmittance is then known at 
each increment of range. As the missile approaches its 
target the atmospheric transmittance increases approximately 
exponentially. If the probability of detection is 
calculated for the seeker system in each incremental 
movement performed by the simulation, the engagement 
scenario probability of detection for the missile would be 
the product of transmittance and original seeker acquisition 
probability generated by its own separate technique. Hence 
if the original Pa were 0.8, and the transmittance were 0.6, 
the total system probability of detection would then be 


0.48. 
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IV. ATMOSPHERIC MODELLING 

The previous chapters presented the principles of 
propogation of infrared radiation throughout the atmosphere 
and the factors which affect the acquisition process. This 
chapter will deal with an empirical approach to modelling 
the atmospheric transmittance using very accurate single 
equations once the conditions of visibility, range and wind 
speed are given. 

Atmospheric transmittance losses for infrared radiation, 
not unlike radio frequency transmission losses, exhibit the 


characteristic 1/(Range)? loss factor as shown below. 


1.0 


Atmospheric 
Transmittance 


ο 3 10 19 20 Range (Km) 


Figure 4-1. Transmittance versus Range. 


An actual curve varies widely in its rate of transmission 
loss with range depending heavily upon the particular 
infrared wavelength in question, the amount and size of 
Particulate matter in the air, the altitude(s), and weather 
conditions as discussed previously. The LOWTRAN program 


CRef. 1] was developed to a large extent from empirical 
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observations and 15 used to predict atmospheric 
transmittance for any given set of conditions. The desired 
goal in the continuing development of the SIREM model was to 
find a simple technique for incorporating the calculation of 
atmospheric transmittance into the incremental range 
movements of a simulated incoming missile and its associated 
infrared seeker head in order to more accurately portray the 


probability of acquisition of the seeker. 


A. )ا‎ 


The LOWTRAN code calculates both atmospheric 
transmittance and radiance throughout the infrared spectrum 
(0.25 to 28.5 micrometers). It uses a single parameter band 
model for molecular absorption, and includes the effects of 
continuum absorption, molecular scattering, and aerosol 
absorption. Refraction and earth curvature are included in 
the calculation for slant paths. 

Five seasonal models and the 1962 U.S. Standard 
Atmosphere are provided as inputs to the LOWTRAN program 
using typical altitude, pressure, temperature, water vapor 
density and ozone density profiles for each. The five other 
atmospheric profiles consist of (1) Tropical (15° N), (2) 
Midlatitude Summer (43° N-July), (3) Midlatitude Winter (45° 
N-January), (4) Subarctic Summer (609 N-July), and (5) 
Subarctic Winter (60° N-January). With these profiles the 


user may select the general atmospheric conditions 
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applicable to the particular simulation. Additionally the 
user may input radiosonde data if desired to make 
calculations for a specific profile. Particular information 
regarding mixing ratios and assumptions made by the LOWTRAN 
program can be obtained from the LOWTRAN Manual CRef. 1] and 
the associated texts listed therein. 

The LOWTRAN code includes three boundary layer aerosol 
types: Rural, Urban, and Maritime. Due to the nature of the 
SIREM/ASM model, only the latter was used since the 
composition and distribution of aerosols of oceanic origin 
are peculiar to the IR/ASM problem. These aerosols are 
largely sea-salt particles in the lower boundary layer 
caused by evaporation of sea-spray droplets. Together with 
a background aerosol of more or less pronounced continental 
characteristics they form a fairly uniform maritime aerosol 
which is representative of the lower 2 to 5 kilometers of 
the atmosphere over the oceans. 

These aerosols should not be confused with the heavy 
direct sea-spray aerosols found in the lower 10 to 20 meters 
above the ocean surface which is strongly dependent on 
recent wind velocities and sea surface turbulence. The 
extinction and absorption coefficients are calculated as a 
function of the relative humidity and based heavily upon the 
refractive index for the given altitude and weather 


conditions. 
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The ability to calculate slant-path transmittances is 
incorporated utilizing a data base of seasonal and latitude 
dependencies of vertical aerosol distributions. This allows 
transmittance calculations of "high-flyer" missile attack 
profiles or aircraft infrared missile engagements by 
changing the input parameters to the LOWTRAN program to 
match those of the simulation requirements. Iterating this 
calculation throughout the slant path range values yields 
the data points necessary to model the reduction in 
transmittance over the entire engagement-range scenario with 
a single equation. The horizontal path used by the SIREM 
model for simulating cruise missiles is assumed to be a 
constant pressure path where earth curvature and refraction 
etfects are negligible, or taken into account through the 
use of the 4/5 earth radius model. An alternative flat 
earth may be assumed for short range situations. 

The LOWTRAN code calculates transmittance as a function 
Of "absorber density" for the path, the pressure, and 
temperature, for the particular wavelength band chosen. It 
utilizes both empirical laboratory data and available 
molecular line constant data in performing its calculations. 
The output for a typical LOWTRAN transmittance run is shown 
in Figure 4-2. The input to the program for spectral range 
is in the form of wavenumbers (inverse centimeters). In 
this example the range is 900 to 1145 cm’ at increments of & 


cm' required by the LOWTRAN program. The lesser number of 
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wavenumbers 15 always the first input into the LOWTRAN 
program, while the higher wavenumber is the second input. 
This corresponds to the high-wavelength to low-wavelength 
output in column 2 of Figure 4-2. The conversion to 
wavenumbers may be made using the following formulas. 
xX - 10927 (WHavenumber s) or 
(Eq. 4.1) 
Wavenumbers = 107*/) 

Here the wavenumbers are in inverse centimeter values as 
used by the LOWTRAN program. Transmittance values for each 
of the contributing absorbers (HO, CO2, Ozone, Ne, Molecular 
Scattering, Aerosol, and Nitric Acid) are then calculated 
for each wavenumber (wavelength in column 2) in row format. 
The total transmittance for each wavenumber is then 
calculated as the product of the transmittances of each of 
the absorbing elements and is listed in column 3. The 
average transmittance listed at the bottom of Figure 4-2 is 
simply the average of the transmittances in column 3, and 
reflects the "integrated" overall transmittance across, in 
this case, the 3.7 to 11.1 micrometer bands for the chosen 
input parameters (model, visibility, range, etc.). 

Since the bandwidth of most IR detectors is limited, 
either by its own material limitations or external filters, 
the input wavenumber range will seldom be as wide as the 
previous example. The 3.5 micrometer band was selected for 


the development of this paper and consists of wavenumber 
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inputs of 2850 (lower) to 2860 (upper) cm ™', and correlates 
to 3.496-3.509 micrometers in 5S cm' intervals. 

In addition to the output shown in Figure 4-2 the 
program generates a large quantity of "header" information 
for each run, which reprints the inputs to the program. The 
unfamiliar user should obtain hardcopy of early runs to see 
what the output consists of. Once comfortable with the 
location of various values of the printout, use of the video 
screen is more appropriate and less wasteful of paper. The 
objective of each run is the "average transmittance” for the 
given input parameters, which will later be used in the 
curve fitting process. 

One can easily see that in order to curve fit the 
transmittance as a function of range with each other 
parameter held constant, enough runs through the LOWTRAN 
program must be conducted at "selected ranges" to obtain a 
smooth curve fit. This means that for every model, range, 
wind condition, etc. the LOWTRAN program must be run in 
order to obtain the appropriate transmittance value. This 
would be a monumental task for the computer, and without 
stripping the PRINT statements from the LOWTRAN code, could 
generate a massive output file. The objective here is to 
show that a minimum number of selected ranges (7 or 8) can 
be used for this process to decrease computer time spent. 
The task becomes much easier if a particular set of 


conditions is in mind, i.e. a Midlatitude Summer model with 
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average wind conditions of 3 m/s and visibility of 50 
kilometers for example. This circumstance could be modelled 
by running a "stacked input" facility of the LOWTRAN program 
containing the ranges desired. A test run can be performed 
to find the approximate maximum range where the 
transmittance value falls below some operator desired 
threshold; e.g. T=0. 0025. 

Curve fitting atmospheric transmittance as a function of 
range requires that all other parameters (1.e. wind, 
visibility, etc.) remain constant. Additionally, enough 
points on the range scale must be taken to provide the 
necessary accuracy across the span of range of interest. 

The missile engagement scenario requires that transmittance 
be a factor in the acquisition process from horizon to 
target, thus the span in range should be from approximately 
<O kilometers to zero. After extensive empirical analysis 
the values in fable 4.1 were selected to meet the above 
criteria. 

Up to this point it was assumed that all parameters 
(including wind and visibility) were constant. The desired 
goal 1s to analyze the individual effects of wind and 
visibility on transmittance and develop an empirical 
equation which generates a good approximation to the 
transmittance value for any wind and visibility condition 
over the scenario range, similar to that shown in Figure 4- 


ave The transmittance takes on a negative exponential 
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characteristic whose decay constant varies as a function of 
wind and visibility. The idea 1s to empirically establish 
the decay constant as a function of both wind and 
visibility, the two principal parameters which must remain 
variable to the model user to realistically model any 
engagement scenario. Intuitively (and verified by LOWTRAN) 
atmospheric transmittance decays much more rapidly if either 
visibility decreases or wind increases. To establish how 
the decay constant changes with wind and/or visibility, 
enough values for both wind and visibility must be analyzed 
to curve fit the decay constant dependence to either. The 
values of wind and visibility in Table 4.1 were empirically 
selected and adequately cover the conditions for most 
scenarios. Then, for each combination of wind and 
visibility, transmittance values over the range span are 
obtained using LOWTRAN. This builds a 3-dimensional table 
which is tha basis for curve fitting and establishing the 
decay constant as a function of three variables; i.e. range, 


visibility, and wind. 


TABLE 4.1. RECOMMENDED CURVE FIT DATA POINTS 


Range in Km Wind m/s Visibility in Km 

552 ο. 1 40 
1.0 oe 5ο 

zu 10-0 20 

4.0 190 15 

8.0 A 10 

120 

18. 0 

21.0 
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The method used for data reduction is an extremely simple 
technique implemented for microcomputer spreadsheets. The 
method employs a scheme for performing Weighted Least 
Squares curve fit techniques to the data, and while 
calculating the coefficients of the curve fit, the goodness 
Of fit of the curve to the data may be viewed graphically by 
the modeller for approval or correction as necessary. Some 
examples of the process will be given later. 

The LOWTRAN program serves as an excellent tool for 
calculating transmittance over a wide variety of parameters, 
however, the program consists of 7000 lines of FORTRAN code. 
Most simulations are already quite large, and the addition 
of this code size as a subroutine is undesireable, 
considering that each increment in range of the missile 
flight path must vector to LOWTRAN as a subroutine in order 
to obtain an updated single transmittance value. 

The alternative method presented here for obtaining the 
transmittance value desired is by utilizing LOWTRAN external 
to the simulation and curve fitting the output data points 
as a function of range to a single equation, given the 
visibility and wind conditions of the particular engagement. 
This method requires extensive use of the LOWTRAN program to 
obtain the data points necessary to perform the curve fit. 


Once completed though, the two coefficients of the 


20 


approximated curve are the only data necessary to write an 


equation of the form: 


T = 9 Exp(—-«FO (Eq. 4.29 


The coefficients of this equation (60 and 4) are all that is 
necessary to allow calculation of transmittance as the 
missile moves incrementally through the stochastic 
simulation. When an increment inrange occurs, a vector to 
the transmittance calculating subroutine (or function) with 
this equation updates the transmittance affect on the 
probability of acquisition. As previously discussed, the 
decay constant («) is actually a function of wind and 
visibility. Later in this chapter accomodations will be 
made for this. 

The SIKEM model specifically calculates through 
simulation, the effectiveness of deployment of Torch/PIP 
Torch type infrared decoys against infrared seeking anti- 
ship missiles (ASMs). It 15 assumed from this point on that 
the incoming missiles utilize a "low altitude” attack 
profile typical of the cruise missile variety which prevents 
early detection by the target ship. Assume for the time 
being that it flies at an altitude of approximately 60 feet 
(20 meters from the LOWTRAN example modelled herein). This 
would constitute a horizontal transmittance path when using 
the LOWTRAN program to calculate transmittance throughout 


the missile-to-target flight path. Some other assumptions 


38 


utilized in the empirical derivations in this paper which 
were used as inputs to the LOWTRAN program are: 
a) The Marine Atmospheric Boundary Layer model was the 
only one utilized. 
b) A land haze factor of 3 was used (default). 
c) The wavelength modelled here was 5.5 micrometers. 
d) Each of the ó atmospheric profiles was calculated. 
e) Default parameters inherent to LOWTRAN were used, 
varying only the visibility conditions, range, and 
wind conditions 


f) There was no precipitation. 


Each of these parameters may be changed if desired to 
meet an alternative modelling requirement. For example, if 
it is known that a "new" IR missile threat has evolved which 
utilizes a "high-flyer" attack profile and a 3.2 micrometer 
optimized detector, the modeller could then choose the 
Slant-path option in the LOWTRAN program selected at the 
necessary wavelength. Knowing only the flight profile, 
detector frequency(ies) and atmospheric conditions would 
then allow full use of the empirical transmittance modelling 
technique presented here without insertion of 9?OOO lines of 
ι ΟΝΤΠΑΝ -—] 

Another advantage to this approach other than its 
simplicity and accuracy is that, if the simulation is 
written in a language other than FORTRAN, this method can 


still be easily adapted due to the short code length 
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required using the standard exponential function. The 
equational form is easily translated as will be seen. 

Let's analyze an example scenario. Suppose we are faced 
with our ship-versus-missile engagement described 
previously, and wish to calculate atmospheric transmittance 
at various ranges in order to approximate its affect on the 
acquisition process. First we must assume some 
meteorological conditions appropriate to our example for 
inputs to the LOWTRAN program. We will select here the 


parameters: 


1) Midlatitude Summer profile. 
2 Visibility is 30 kilometers. 
5) A maritime environment. 


τν 


4) A land haze factor of 3 (may vary from 1 to 10, with 
10 being bad haze as around Los Angeles). 


5) Wind is virtually nill. (0.1 m/s) and has been so 
for the past 24 hours. 


6) The range from the seeker to target is 4 kilometers. 
7) Assume a 3.3 micrometer band (28530-2860 cm’) 


8) The path 1s horizontal at 20 meters altitude. 


We may now insert out data into the LOWTRAN program in the 
format listed in the LOWTRAN Manual (Ref. 11. 

In doing so we would find that at our 4 kilometer range, 
the transmittance for this wavelength would be 0.63516 (See 
Appendix A: fable A.2.a). The value may differ only 


Slightly between computers. This means that the received 
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energy density at a range of 4 kilometers is only 65% as 
high as that for a perfectly transmitting medium (vacuum). 
We can now duplicate our efforts for the 4 kilometer range, 
but now with a wind velocity of 10 meters per second, 
keeping visibility the same. Doing this we obtain a 
transmittance value of 0.4537. In the diagram below we see 
that increasing the wind, or decreasing visibility, has a 
dramatic effect on the decay constant associated with the 


transmittance. 


1.0 
Atmospheric Wind 
Transmittance 


O m/s 


Wind 


10 m/s 


Range (km) 


Figure 4-2. Wind Effects on Transmittance. 


Note the approximately exponential fall off with range. 


This is the curve we wish to approximate as: 


T(R,V,W) = G-EXP (X (CD) - F CHO - R- ÓD (Eq. 4. 9) 


from the data points obtained from the LOWTRAN program. 
Here T, R, V, and W represent transmittance, range, 
visibility, and wind respectively. There is correlation 
between wind and visibility, and the effects are accounted 


for in the factor 4. 
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In order to model any particular atmospheric profiles, 
the data points for the 3-dimensional table may be obtained 
by adapting the following overlay algorithm as the main 


program, using LOWTRAN as a subroutine. 


FOR Atmospheric Profile —- 1 to & * See Note 1 
FOR Wind Speed = 0 to 20 Step & * See Note 2 
Wind = 0.01 -* FLOAT (Wind Speed) 
FOR I = 1 to Nr of Range Values 
READ Range Value 
OPEN Input File * See Note 3 
WRITE Input File 
CALL LOWTRAN 
WRITE Qutput 
NEXT I 
NEXT Visibility 
NEXT Wind Speed 


NEXT Atmospheric Profile 


Note 1: Frior to entering the looping structure the user 
must assign values to the other parameters used by the 
LOWTRAN program (e.g. wavenumbers, haze factor, etc.). 

Note 2: The actual wind speed may not be zero. If so 
LOWTRAN defaults to 4.1 m/s. Hence the following statement. 
Note =: The input file must be formatted as specified in 


the LOWTRAN Manual tRef. 11]. 
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Again it should be emphasized that the source code of 
the LOWTRAN program should be commented out to prevent 
printing of the massive administrative data. Only the 
"Average Transmittance" over the sufficiently narrow band is 
desired. The output may be directed to a printer or file. 
The latter method is preferred if input to the curve fitting 
process is desired ina direct fashion. In order to keep 
track of which data value is associated with which range, 
visibility, and wind value, it may be wise to write those 
associated values in the same output row. 

Continuing with our example we will assume for now that 
only one value for wind, and likewise visibility are 
calculated. The transmittance would then become only a 
function of range since all other parameters remain 
constant. The desired goal is to select range values which 
best facilitate the curve-fitting process. Using the range 
values from Table 4.1 we can now perform an example of the 
curve fitting process, while presenting the spreadsheet 


technique used for doing Weighted Least Squares. 
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2 a ত 


Making the assumption that all parameters remain 
constant as before, with only range changing, a curve fit of 
the following LOWTRAN data is performed for the Midlatitude 
Summer profile with wind at 5 meters per second, and 


visibility of 30 kilometers. 


TABLE 4.2a. LOWTRAN DATA FOR CURVE FITTING 


RANGE TRANSMITTANCE WEIGHTING FACTOR 
OI 0-9272 8 
1.0 ο. 7824 1 
2.0 Oo. 6588 2 
4.0 0.4828 21 
8.0 Oe 27 25 & 
12.0 0.1588 1 
18.0 0.0727 1 
2150 O. 0498 Z 


The Weighted Least Squares procedure is conducted in 
two parts. To effect the weighting, implementing multiple 
occurrence of the same data points is shown in table 4.3. 
Following in figure 4.4a and 4.4b is the actual Least 
Squares algorithm implemented in spreadsheet format. The 
output desired is the values for the A-Coefficient (9) and 
B-Coefficient (ας) for equation 4.2. In this example the 
first coefficient 0-0.875824, and the second, «=—-0. 140899. 
Using equation 4.2 yields the table which follows at the top 


of the following page. 
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TABLE 4.2b. DATA POINTS EXTRACTED FROM FITTED CURVE 


RANGE TRANSMITTANCE 
0.2 0:8515 
1.0 0.7607 
2.0 0. 6607 
4.0 0.4985 
8.0 0.2837 

120 0.16195 

18.0 0.0693 

21.0 0.0454 


Comparing table 4.2a with 4.2b (i.e. the LOWTRAN 
data versus the exponential Curve-Fitted data) we notice 
comparable transmittance values, especially at ranges 
greater than one kilometer. Due to the nature of the 
engagement scenario, we may assume transmittance values at 
ranges less than one kilometer or so to be equal to 1.0, 
since its effect on acquisition at these ranges is minimal. 
Above 1 kilometer the error is generally one to two percent 
at most, and hence, very useable for "probability" 
modelling. Within the following tables the Least Squares 
process in implemented into spreadsheet format. Although in 
this case Lotus Development Corporation's 1-2-3 (TM) was 
used, the technique is similar if using commercially 
available packages. The data is arranged in row-versus- 
column format with the Least Squares algorithm implemented 


in Tables 4.4a and 4.4b. 
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—_ ce em eee ee ee eee a سس‎ 
کک پت س‎ a ee রা ee ee ee ee ee ee ee ee ee ee مس‎ পানা 


I Model: Midlat Sum 
2: Visibility: 30 
3: Varied param: Visibility 


=>: Range (km) Transmittance 


: 029 ο στο 
z 0.2 0.9272 
8: 0.2 0.9272 
: 0.2 0.9272 
10: 0.2 0:9272 
11: 0. 2 0:9272 
12: ο.» 07722 TABLE 4.3. 
15: 0.2 009272 
14: 1 O. 7824 TRANSMITTANCE AND RANGE 
153 2 O. 6588 DATA POINTS USED 
16: 2 O. 6588 
17: 3 0.4828 
18: 4 0. 4828 
12: x 0.4828 
20: E Ο. 4828 
21: E 0. 4828 NOTE: 
22: c O. 4828 
25: + 0. 4828 Using the same data point 
24: 4 0.4828 multiple times at any point 
293 4 O. 4828 minimizes error about that 
26: 4 O. 4828 point in the curve fit 
27: d 0.4828 process. This emulates the 
28: E 0.4828 weighting matrix commonly 
29: > 0.4828 used in performing Wieghted 
S0: 4 O. 4828 Least Squares curve fitting. 
351: 4 Q. 4828 
σος 4 0. 4828 
35: + 0.4828 
354: 4 O. 4828 
d 4 ο. 4828 
56: 4 ο. 4828 
57: 4 0.4828 
38: 8 0.27235 
Og 8 22272 
40: 8 Of 2725 
41: 8 0.2725 
405 8 (০2712 
43: 8 ج272 ے0‎ 
44: 12 0.1588 
45: 18 9-0727 
46: 21 0.0498 
47: 21 0.0498 
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TABLE 4. 4a. 


CALCULATIONS USED IN LEAST SQUARED FITTING 


: C D E £ F 
Se ED am Tn سے جج ۹ ہے‎ Ea سے سے‎ তু এ SS SS SS SS Se SS SS SS SE 2 ےج ے‎ EEE == == 
1: SUM of R, This section finds the coefficients 
: à and B for the equation: 
: 210.6 T - C*Exp(«X*R)..... at Visibility =30 km 
5 Πι2 LntT,) CLn (T,) J2 Rı#Ln (T, ) 
: 0.04 —0. 0796 0.9057 -0.0151 
p 0.04 —0. 0756 0.0057 -0.0151 
: 0.04 —-0.0756 0.0057 —O. 0151 
: 0.04 —0. 0756 0.0057 -0.0151 
10: 0.04 —0.07956 0.0057 =0. 0151 
IMs 0.04 —0. 0756 0.0057 —0.0151 
12. 0.04 -Ο.Ο756 0.0057 —O. 0151 
es ο. 08 -0.0756 0.0057 -ο. ΟΙ 51 
14: 1 -Ο. 2454 0.0602 —0.2454 
15: 4 -0.4173 ο. 1742 -ο. 8.547 
16: 4 —-O. 4173 0.1742 —0.8:547 
07 16 -ο. 7282 0.5302 -2.9126 
18: 1০ z0:7282 0.5302 =2 9126 
po 16 -Ο. 7282 0.5302 —2.9126 
20: 16 -Ο. 7282 0.5302 -2.7126 
21: 16 το. δες 0.5302 -72.7126 
Z2: 16 -Ό. 7282 0.5302 ০৪5 
Bos 16 -ο 7282 0.5302 —2.9120 
24: 16 πο 282 0.5302 -:229126 
2s 16 —0.7282 0.5302 -2.9126 
26: 16 205721927 0.5302 2590726 
27: 16 -ο. 7282 0.5302 EDIT 
28: 16 -ο. 7282 0.5302 2209126 
29: 16 -Ο. 7282 0.5302 2ج‎ 60126 
30: 16 -Ο. 7282 0.5302 2 αν σα 
3i: 16 -Ο. 7282 0.5302 —2.9126 
oz: 16 -Ο. 7282 0.5302 চির 5 
E. 16 —0.7282 0.5302 2-9126 
34: 16 πω 7282 0.5302 —229]26 
Sos 16 —O. 7282 0.5302 —2.9126 
36: 16 -ο. 7282 0.5302 —2.9126 
57٤ 16 -Ο. 7282 9D. 5302 72.9126 
38: 54 —1. 3009 1.0922 —10. 4068 
E 54 -۱ 3009 126922 —10.4068 
40: 64 -1.3009 1.6922 -10. 4068 
41: 54 —1 . 3099 1.6922 -10. 4068 
42: 54 —-1.3009 1.6922 -10 68 
43: 54 21235009 1.6922 —10.4068 
44: 144 —1.8401 3.5860 —22.0813 
45: 324 —2.6214 6.8718 —47.1855 
46: 441 =2. 9997 8.9984 —62.9945 
47: 441 "+8 9 8.9984 —62.9945 
48: Sum (R, )? Sum LntT,? Sum Squares Sum of Products 
49: 2079. 35200 —35.2421 49.9966 —320.8972 
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TABLE 4.4b. SPREADSHEET LEAST SQUARES CALCULATIONS 


: G : H : Formula in cell at left: 
========================== Note: Here n=42 
Scratchpad : 
—0.1409 2 (42*F49—-D49*C.5) /(42*C49—-C35^2) 
=0. 1526 > (D49—-G2*C3)/42 
Μα δι ος 7 62#(F49-C3*D49/42) 
20.4251 > E49-D49^2/42 


A-Coefficient 0.875824 > Exp (63) 
B-Coefficient -0.140899 » G2 
R-Squar ed 0.9946 > 64/65 


৭০ 00 এ 0 01 ν»» 


This satisfies the Least Squares process. Cells G2 through 


GS represent respectively: 


51) ΓπΣίΒι]πίτι))} - El1In(T))ER,J/EnER;? - (ΣΒΙ}:1 
G2) [5:17 701) - )61( ) 0٥71ہ‎ 

(Eq. 4.4-4.7) 
G3) 62 LER. 1n(T,) - ER,Eln(T.) /n3 


G4) τι) πο στο το επ 
Cell C3 is the sum of the range values of column "A". 
C3) ER, (Eq. 4.8) 


The R-Squared value represents the "Goodness of Fit" 
of the function to the original data when using the 
approximating equation. In this case the fit is a very one, 
but in cases of either extremely good or poor visibility 
(i.e. *40 or <3 kilometers) errors between LOWTRAN 
predictions and empirical results grow, and the fitted curve 
will have regions where error may be on the order of five to 


Six percent. For modelling in the probabilistic sense this 
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is still adequate, since the region of error is transient 
and errors of this magnitude have only small effects on the 
acquisition probability. Figure 4-3 is a plot of the above 
example using both LOWTRAN-obtained data points, and the 
corresponding exponential curve fit data points obtained by 
using the spreadsheet technique presented. We may now 
extend our process to calculate the atmospheric 
transmittance for any of the 6 LOWTRAN profiles for any wind 


or visibility conditions. 


Figure 4—.5 


Comparison — LOWTRAN vs Exponential Fit 


Atmospheris Transmittance 





ο 2 4 6 8 10 12 14 16 18 20 
Range in Kilometers 
0 LOWTRAN Data + Curve Fit Data 
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a. Extending the Curve Fit 


The 3-dimensional table (A.2.a) of transmittance 
Values has been obtained for, in this case, the Midlatitude 
Summer Profile. Thus far we have completed the curve- 
fitting process for one particular wind and visibility 
condition (1.e. visibility of 30 Km, wind speed 5 m/s as 
calculated in Table 4.4). The A and B coefficients (60 and 
a), which were derived as a result of Least-Squares fitting 
an exponential function, were used in equation 4.2 to 
emulate extinction over range for this wind/visibility 
combination to obtain Figure 4-3. There the accuracy of the 
exponential fit versus the original LOWTRAN data is visually 
apparent. Tables A.l.a through A.6.b in Appendix A list the 
Values of transmittance obtained from the LOWTRAN program 
and from the curve fitting process for each of the six 
atmospheric profiles available. The curve-fit coefficients 
for the other five atmospheric profiles are contained in 
Tables A.7 and A.8 in Appendix A. 

We now wish to extend this procedure to be able 
to emulate extinction over range for any wind/visibility 
combination. To do so we must first build a table of the A 
and B coefficients and examine how their values change from 
one wind/visibility combination to another. Using the 


combinations from Table 4.1 the coefficients were obtained 
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for the entire atmospheric profile of our example. These 
are contained in Table 4.5. 

By examining Table 4.95 we see that the A 
coefficients hardly change throughout entire range of wind 
and visibility values. 


This coefficient may be assigned a 


constant value for the remainder of the extinction curve-fit 


process. This is true when modelling each of the 
atmospheric profiles individually; however, the A- 
coefficient does change from one profile to another. This 
can be seen in Table A.7 in Appendix A. 
TABLE 4.5. EXPONENTIAL CURVE FIT COEFFICIENTS 
FOR "MIDLATITUDE SUMMER" PROF ILE 
A—Coefficient Value 
Vis=40 Vis=30 Vis=20 Vis=15 Vis=10 
eee he EEE ee ee ৯১৯৯ M 
0 0.875788 0.875808 0.875734 40.875854 0.875812 
9 0.875747 9.875824 0.875699 0.875477 0.877011 
10 0.875820 0.875716 30.875650 0.874656 0.874052 
z 0.875753 0.875685 0.875366 0.876011 0.872420 
20 0.875777 0.875885 0.875886 0.877431 0.859462 
B-Coefficient Value 
Vis=40 να 5 ο) Vis=20 Vis=15 Vis=10 
০ তা 22 επ 
Q —0.064275 —0.065945 —0.06931.3 —0.072711 —0.079494 
J —0.118579 —-0.140899 -0. 185840 -0.231112 -0.322254 
10 =O οσους -0. 1565973 -0.207757 -0.263695 -ü. 57/1488 
Io —0.1.2919 —-0.1606635 —0.216477 —0.273163 —0..:84626 
20 —0.1:54594 0.162714 -0.219772 -ü. 27/7832 —0. 386820 


The B-coefficients do vary throughout Table 4.5. 


Through further examination, the change in its value is non- 


linear with wind and visibility, and as such, simple linear 


1 


interpolation of this table is not accurate. Nor is this 
technique desirable because it requires many more B- 
coefficient data points than those contained in Table 4.5 to 
be stored in memory to be even moderately accurate. The 
process of extending equation 4.2 to the more general form 
of equation 4.3 becomes a matter of finding out how the B- 
coefficients in Table 4.5 change as a function of wind and 
visibility (i.e. find the exponential decay constant of 
extinction in terms of its dependence on wind and 
visibility), and fit functions which will accomodate the 2- 
dimensional decay constant variation. 

Thus the curve fit process is continued, only 
now it must be performed on the B-coefficients which vary 
with both wind and visibility. This will leave ‘Range’ as 
the only variable to input to equation 4.5 for calculating 
transmittance. From table 4.5 it can be seen that the B- 
coefficient increases negatively (decreases) as either 1) 
wind increases, or 2) visibility decreases, and is therefore 
a function of wind and visibility. Equation 4.3 stipulates 
that the decay constant for extinction over range consists 


of three parts. These are: 


8) αν) — Due to changes in visibility, 
b) f(W) — Due to changes in wind, and 
C) ó - Any correlation between the first two. 
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By also observing that the B-coefficient is 
least negative where wind is minimum and visibility is 


maximum, we can now use the ratios of any other B- 


coefficients relative to the first and find their 


multiplicative increase over the original (baseline - 


.------ 


maximum visibility, minimum wind) B-coefficient due to wind 


and visibility. This done by dividing each value in the B- 


coefficient table by the ‘baseline’ value in the upper left 


corner of the table as is done for Table 4.6. Table 4.6 


depicts how much greater (negatively) the decay constant is 
for any particular wind and visibility combination than it 
would be for its best transmittance (i.e. wind zero, 
visibility maximum). The non-linearity of the B-coefficient 
with wind or visibility is much more evident in Table 4.6. 


We shall see in a few moments why this table is necessary. 


TABLE 4.46. RATIO OF BiVis,Wind) TO B(V1 Snax, Windarn) 


Vis Vis Vis Vis Vis 

4ο 3 20 15 10 
EL ul le Lu uuu LL LL Li σσ سے‎ 
Q 1.0 1.02598 1.07858 2 1.23678 
D 1.84487 2.19215 2.909153 35.509567 οι οι ον 
10 2.02035 2.435322 5.26657 4.10257 J. 77966 
15 2.06797 2.49962 AX. 56801 4.24991 3.98407 
z 2.09092 292195 3.41924 435225 6.01820 


The following steps are the procedures used to 


obtain the three functions of the exponential decay constant 


of equation 4.5. 


functions «(V), 


f(W), 


Individually these steps produce the 
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and ὁ which will allow the input of 


any wind or visibility value from 0-20 m/s or 40-5 Km 
respectively, and generate the approximated B-coefficient of 
extinction for that condition. Keep in mind that this 
example only covers the ‘Midlatitude Summer’ profile. The 
other profiles are approximated in the same fashion, with 


the results included in Table 4.7. 


(1) Step One. As visibility decreases the 
B-coefficient becomes more negative, causing transmittance 
to decrease at a faster rate (more extinction). Testing of 
various functions was conducted to generate the B- 
coefficient as a function of visibility. It was assumed to 
Start with that wind is zero so that the curve fit process 


uses the B-coefficients of the first row on Table 4.5. 


Graphically the data points appear as in Figure 4-4. 


Figure 4—4 


B-—Coefficient Curve Fit to Eq. 4.9 
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Care must be taken in selecting a curve-fit 
method to approximate any of the three functions «(V), f(W), 
or 46 since the B-coefficient data points are negative. This 
eliminates exponential or logarithmic function Least-Squares 
techniques since the natural logarithms of the data points 
are necessary, and are undefined. 

To approximate «(V)», a cubic function was 


selected of the form: 


Row E«í(Vis)]J: A, + Az2-V + As-V*% + Age V® (Eq. 4.9). 


This eliminated the negative data point problem and produced 
highly accurate results for the few data points on hand. 

The values found for the coefficients of the cubic are 
contained in Table 4.7. These values were used to generate 
the fitted curve in Figure 4-4, while the original B- 
coefficients were plotted as the single point values. 
Remember: fitting of this function is performed only on the 


first row of B-coefficients in Table 4.5. 


(2) Step Iwo is to establish a function which 
reflects the change in the decay constant with wind: i.e. 
f(W). Graphically the data points plot as shown in Figure 
4—5. At first the function appears as an exponential decay 
which would be easily modelled. But, since the data points 
are negative this method is not easily performed. This is 


where the concept of Table 4.6 is useful. Given that the BR- 
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Figure 4—5 
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Actual B-coefficient change with increasing wind. 


Figure 4—6 
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B-coefficient translated for curve-fit to Eq. 4.]0. 
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coefficient for any visibility condition (with wind zero) 
can be generated by completing step one, then f(W) becomes 
merely a multiplier to that value generated. 

To find out what to multiply a(V) by to get 
the decay constant for wind conditions which are not zero, 
we fit a function to the first column of Table 4.6. Now, 
all the data points are positive, and any function is fair 
game. Graphically, the newly translated function appears as 
in Figure 4-6. After unsuccessful (inaccurate) attempts to 
fit a hyperbolic tangent (tanh) function to this curve, a 


much more accurate one was discovered: 


Column Cf (Wind) J: 1+€D,- (1-Exp (Dz-Wind))] (Eq. 4.10) 


(3) Step Three. Figure 4-6 shows a family of 
five curves which plot each of the five columns of Table 
4.6. The curve with the least maximum represents the first 
column in Table 4.6. The others in increasing fashion 
represent decreasing visibility conditions. The correlation 
factor ó is necessary to increase the extin-cion when wind 
and visibility are not zero and maximum respectively. In 
this case 6 acts as a multiplier also, driving the decay 
constant more negative when its value becomes greater than 
one. Fitting a function to points vertically spaced on 
Figur 4—6 is the key to finding the correlation coefficient 


Bs A function which works extremely well in this case is: 


S7 


c (Eq. 4.11)‏ تہ 


The same effect may be obtained by fitting 
equation 4.11 to the data points of one of the lower rows of 
Table 4.6 if they have been normalized a second time in the 
horizontal direction (i.e. divide each coefficient in one 
lower row of Table 4.6 by the coefficient in the first 
column of that row). 

This procedure completes the curve fitting 
process for each atmospheric profile, leaving eight constant 
coefficients which, when incorporated into equation 4.3, 
yield an admirable approximation of transmittance values at 
range values between one and twenty kilometers throughout 
the entire atmospheric profile. This procedure was used to 
establish the constants for each of the six LOWTRAN 


profiles. These constants are included here for reference. 


TABLE 4.7 


CURVE FIT COEFFICIENTS FOR ALL SIX LOWTRAN PROFILES 


TROPICAL MIDLATITUDE MIDLATITUDE 
PROFILE SUMMER WINTER 
ων ; -0.139154 =0. 12475992 το 087365 
Az ¦ 6.2887E-035 62796106 7 6. 31735317E 05 
As τ —2.23554E-04 —2.207060E-04 —2.247950E-04 
^ : 2.6091E-O6 2.60194B82E-O06 2.6:777/25E-06 
C, +» 15.741035 15.182516 29.: 876 
C2 : -Ό. 7168046 —0. 743405 -0.925837 
D: 0.895 1.10 2.46 
آبە+ و0‎ IA τος AS 
0 ; ο. 88Ο 0-875 0 5 | 
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A, 
Az 
As 


C2 
D, 
Dz 


SUBARCTIC 
SUMMER 


EO T0973: 
Ó. 25456-03 
=2. 22286-04 
2e IISSE—-O6 
pou qz] 
—0.80492 
1.440 
πο 0 

ο. 


TABLE 4.7 


SUBARCTIC 
WINTER 


—0.07970 
6.4225E-0S 
-2.28097E-04 
2.60725E—O6 
57.730 
—0.98561 
5.980 

SO: 22S 
0.9636 


(Continued) 


STANDARD 
ATMOSPHERE 


—0. 091037 
cia Qc -0.5 

—2.0021E-04 
2. o5407E-06 

24.130 

—0. 86706 
IBI 

—0. 354 
0.9208 


Note here that eight of the nine constants 


in each atmospheric profile eventually constitute «, the "Ε΄ 


or extinction coefficient. 
and explain the computer coding process. 


contained in appendi: 


The next step is to implement 


This procedure is 


We can now turn our attention to 


the SIREM model and examine the acquisition process. 
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V. THE STOCHASTIC INFRARED ENGAGEMENT MODEL (SIREM) 

As discussed in chapter one, the most inexpensive means 
of developing equipment characteristics and tactics for 
their use is through simulation. The use of computers to 
analyze effectiveness of systems or subsystems is becoming 
the predominant development method for this reason. The 
SIREM model was designed to investigate the effectiveness of 
deploying decoys which emit an infrared signature against 
infrared seeking anti-ship missiles (ASMs). 

By examining through simulation such parameters as burn 
time, deployment altitude, deployment angle, decoy drift and 
descent characteristics if airborne, and distance from the 
ship, optimum design requirements and deployment tactics can 
be developed. As with any software model, to maintain its 
Validity it must be periodically reviewed and upgraded to 
remove any synthetic qualities and replace them with 
accurate representations. 

There are ongoing efforts to incorporate atmospheric 
losses, discriminative seeker capabilities, and last but not 
least the multi-ship multi-missile analysis capability. 
Coincident with the work being done on the SIREM model, 
extensive investigation is being conducted regarding optimum 


deployment characteristics for chaff, repeater decoys, anda 


60 


variety of other offboard countermeasures in order to 
integrate the use of offboard devices into a multi-band 
defensive system. To fully understand the SIKEM model 
requires extensive programming experience and/or simulation 
expertise. Likewise an operational knowledge of available 
devices is advisable; however, a general description of the 


program is in order. 


A. THE MODEL 


The SIREM model calculates probabilities associated with 
distraction of an IR homing missile by a pattern of decoys 
launched by the attacked ship. Rather then the Monte Carlo 
method of simulation, a stochastic approach was utilized in 
order to significantly reduce computation time. The program 
uses realistic input parameters pertinent to the scenario 
under study, and evaluates statistical information about the 
distraction effectiveness which may be output in graphical 
Or numeric form. 

The user may easily modify any of the input parameters 


under the following general headings: 


1. Ship characteristics 
2. Missile characteristics 
5. Decoy characteristics 


4. Scenario and run characteristics 
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Although the original version transported from Naval 


Research Laboratory was designed to operate under the NLX 


NAMELIST facility of the Prime (TM) computer there, 


significant changes were required to adapt the input/output 


to non-NAMELIST type since it was not available at Naval 


Postgraduate School; and to modify the graphics subroutines 


and 


the 


hardware to operate on the Digital VAX 11/780 (TM) and 
DI-3000 (TM) graphics software available there. 


Output from the program is available in the form of 


graphs or printed output. A summary of the available graphs 


15 1 


isted here: 
Frobability of distraction conditioned on initial 
missile-to-ship distance. 


Frobability of distraction conditioned on initial 
azimuthal approach angle. 


Fattern of ship lock-on range as a function of angle. 
Frobability of distraction conditioned on initial 
azimuthal approach angle and a single initial missile 


to ship distance. 


Angular regions where probability of distraction 
exceeds a threshold. 


Locations where threshold crossings of plot S occur. 
Probability of distraction conditioned on initial 
missile-to-ship distance and sector of initial 


azimuthal approach angle. 


Ship orientation, decoy deployment, and wind 
direction. 


Effects of atmospheric refraction on IR at various 
incidence angles. 
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Being stochastic by design, certain parameters may be 
considered random variables for initialization. The initial 
position and orientation of the missile relative to the ship 
is a good case. It is assumed for startup that the missile 
is pointed radially inward toward the ship, and hence the 
effectiveness computations would yield a lower bound for 
effectiveness of decoy distraction. This is the basis for 
the 'worst case' situation described in earlier chapters. 

The initial range and azimuth approach angle of the 
missile are discrete random variables. The initial range 
may be changed to appropriate kilometer values; while the 
initial range may also be modified, it is designed to run 
multiple missiles at the ship every S60/Ninera degrees. 

The seeker characteristics and the random variables 
assigned to emulate it can change from missile to missile. 
In this sense, randomness must be assigned to accomodate 
scanning directions, boresight error, field of view, and 
probability of acquisition on multiple targets (i.e. ship 
plus any decoys). The model performs these adequately, 
using significant but valid assumptions regarding the 
process. 

The maximum initial range of the missile to the target 
is R at the time of launch of the first decoy. This is 
considered time t=0. It then takes Ro/Vn seconds for the 
missile to reach the ship at the origin, with V, being the 


missile velocity. The sample time T, interval determines 
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the maximum number of time samples taken Ns. T, is set 
small enough to obtain an accurate representation of the 
scenario. Since the missile is always aimed at ship before 
acquisition, a lower bound on distraction results, and if no 
decoy is acquired then the probability of distraction is 
zero for that given initial range and azimuth. 

The choice of time interval T, is dependent on the 
missile velocity, target velocity, and the scanner sweep 
time Ts. Since the area of the scanner field of view on the 
Sea surface is large relative to the distance either the 
missile or ship travels during the sweep time, Ts may equal 
"Eus As different conditions arise a necessity for time 
scaling may result. Developments have begun to incorporate 
different scanner parameters. 

The projection of field of view is incorporated to allow 
both curved earth and flat earth calculation. The ship and 
decoys are treated as point targets for infrared radiator 
purposes. The "hot spot" of the deployed decoy is placed at 
a known distance above the sea surface. Calculations are 
made to determine if the hot spot is within the field of 
view of the scanner. The next stage of development requires 
incorporating radiance and transmittance calculation, and 
detector sensitivity parameters to determine if there is 
sufficient radiance from the ship or decoy to be detected by 
the seeker. Much of this development was covered in earlier 


chapters. 
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The model can also calculate any shadowing of decoys by 
the ship. This means that if the aspect of the ship comes 
between the missile and decoy, effectiveness of that decoy 
will inherently drop to zero. This can be done using either 
the flat or 4/3 earth model. In addition to shadowing by 
the ship, if excessive range exists, shadowing by the earth 
may occur. 

Once acquisition of a decoy occurs, the missile flies 
according to the aerodynamic parameters of the control 
surfaces and systems for whatever missile is being modelled. 
If a decoy is not in the field of view and the ship is 
acquired, the missile is assumed to hit the ship. If locked 
onto a decoy, the missile flies at the decoy until it passes 
it. At that point break-lock occurs. This condition may 
also occur if the control system forces the missile to 
maneuver so that the decoy is no longer in the field of 
view. find lastly, if the decoy burns out break-lock occurs. 

When break-lock occurs the scanner must re-enter the 
acquisition mode. Since the missile usually makes drastic 
maneuvers near break-lock, the re-acquisition process is 
somewhat random, but if the ship is in the field of view it 
may be assumed to hit the ship. Missile flight continues at 
break-lock until all targets are passed. 

The flexibility of the program is excellent. Virtually 
any parameter is changeable either at input or at source 


code level. As many as five decoys may be deployed by the 
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ship (more if certain array dimensions are changed) in any 
combination of deployment angles, times, altitudes, and 
active lifetime. The ship may change size and speed. The 
wind is a factor on decoy drift and missile flight. And the 
missile may change speed, altitude of approach, and control 
system or seeker parameters. The variations are endless, 
and analysis of the many combinations pertinent to real- 
world scenarios may be undertaken to make maximum effective 
use of offboard countermeasures. Additionally, by further 
extending not only SIREM but other similar models as well, 
evaluation of descrimination techniques and optimization of 
decoy signatures become possible avenues of further 
development. With this in mind we turn our attention to the 


acquisition process and observed phenomena. 
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VI. IR SEEKER SYSTEM AND COUNTERMEASURES 

Since the principal interest here lies in passive 
seekers, it 15 convenient to separate these into three broad 
categories: (1) Scanning systems such as search, track, and 
imagers, (2) Staring systems which use non-mechanical means 
to obtain spatial scene viewing, and (35) Measurement systems 
such as radiometers, spectrometers and interferometers. Of 
interest here are the first two categories. Scanning 
detection schemes are most common since staring systems 
usually require sophisticated multi-dimensional array 
detectors and massive numbers of amplifiers to support the 


individual elements. 


A. SCANNING AND STARING SYSTEMS 


Scanning systems sample the radiant intensity 
distribution in their field of view using various 
horizontal/vertical techniques. The output from the 
detector in this case will be a linear analog of the 
distribution or a simple indication of the presence of a 
target and its location. The object plane is usually 
assumed at infinity and the detector responds to radiant 
levels within the instantaneous field of view of its 


element (s). This gradient in radiance is translated to 
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electrical signals which are then filtered, amplified, and 
referenced to a baseline (missile axis) for inputs to the 
various azimuth and elevation control systems. For imaging 
systems the signal is translated from the object to its 
corresponding image plane by painting out the radiance 
distribution using optical deflection onto a film, or beam 
deflection to a cathode ray tube. Using internal processing 
techniques in missile seekers, the translation is used as 
input to the missile control surface system to center the 
highest radiant intensity in the center of the total field 
Of view. The use of either spiral or raster techniques to 
accomplish slow mechanical scanning is most common. In some 
cases the relatively high speed motion at altitude provides 
the vertical trace capability for the system, while a 
spinning mirror provides the horizontal trace. Regardless 
of the scanning technique employed, major concern has 
recently been given to the incorporation of discrimination 
processing into newer seeker systems. If (or, more 
appropriately, when) this occurs, the effectiveness of using 
offboard decoys could be significantly reduced. Therefore, 
an examination of possible discrimination methods is in 
order. 

In dealing with imaging systems there are two basic 
approaches to describe the sensitivity of the detector. Ihe 
first approach is based on the detector figure of merit (D*) 


described earlier. It regards the detector as a noise 
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source of fixed level which is determined from the 
measurements of fixed conditions similar to those 
encountered by the detector when installed in its normal 
operating configuration inside the seeker. 

The other approach which is suitable only for BLIP 
(Background limited mode) detectors regards a random arrival 
rate of photons from the background as the ultimate source 
of detector noise. In such a case the detector noise power 
density is described using the laws of photon rate and 
detector quantum efficiency. When the number of incoming 
(target) photons is small in comparison to the number of 
background photons the approaches are equivalent, since BLIF 
D* 15 determined by background photon rate. If the number 
of target photons is larger than the average number of 
background photons though, the noise due to the random 
arrival time of target photons dominates and D* then is 
meaningless. This case is the unusual one. 

Image forming scanners are not usually designed to 
respond to an absolute radiance level, i.e. not for a dc 
response. They are designed to follow changes in radiance 
as the field of view is scanned. This includes gradual 
changes in the radiance level corresponding to low 
frequencies. DC restoration can be used to control the 
intensity levels at the output. Without this form of 


automatic gain control (AGC) periods of signal blackout or 
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saturation would occur due to high 1/f noise and amplifier 
drift at low frequencies. 

The signal of interest from the scene is the variation 
in irradiance H caused by the variation in temperature or 


emissivity of objects in the scene and may be expressed as: 


00 D 
AH —|w| {[ATfoW dX + jus Ws dX in Watts (Eq. "تک‎ 
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T = temperature in degrees Kelvin 
€, = emissivity (wavelength dependent) 
W = spectral radiant emittance in W/cm? 


If most noise in the system is caused by the detector, 
it may be expressed in terms of "noise equivalent flux 
density” (NEFD). This is the minimum signal irradiance 
capable of producing a peak signal to rms noise ratio (Sp/Nr) 
of one. 

The ultimate goal is to find the actual signal-to-noise 


ratio (S/N) generated by the scanner. This 1S given as: 


S = Variation of irradiance in scene (AH) in Watts (Eq. 6.2) 


ne n ——‏ پیم سے 
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N NEFD 


The noise equivalent power (NEP) in watts of the system 
can be obtained from the figure of merit of the detector 


(D*) and is expressed as: 
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NEP = VAg- 4f in Watts (Eq. 6.3) 


p“ 
where Aa = area of the detector in ۶ 

&f = system noise bandwidth 

D* = detector figure of merit in cm-Hz'/^/watt 


The signal power on the detector is equal to the optical 
collector area times the signal flux density. Therefore the 


system noise equivalent flux density is: 


NEFD = 4 VAg- Af in Watts/cm* (Eq. 6.4) 
8021(۳ 
where D, — diameter of the aperture of the optical system 


eti 
Il 


o optical efficiency 
The difference in temperature and emissivities (AT and 
Δε) of objects in the scene contribute to the differential 
radiant emittance (AH), but due to the small contribution of 
Δε for non-BLIF detector operation, its term in equation 6.1 
is commonly neglected. By substituting Wy, for the first 
integral term the abbreviated value for radiant emittance 


becomes then: 


AH = W-ATW, in Watts/cm? (Eq: 6-3) 


The instantaneous field of view W is defined by the 


detector area ìa and the focal length fl as: 


w= Ag in steradians (Eq. 656) 


آ7 


By substituting equation 6.6 into equation 6.4 the noise. 
equivalent flux density is reestablished as: 
NEFD = AF A/w-Af in Watts/cm? (Eq. 6.7) 
το ρε, 
where F is the f-number of the optics. 


The signal-to-noise ratio generated by the scanning 


System is then (from equation 6.2): 


S = ATW D,D*c,Vw (Eq. 6.8) 
N AF VAf 


Some careful considerations are in order here since the 
AH actually seen by the sensor is modified by the 
transmittance of the optical components and the atmosphere 
as discussed earlier. If the bandwidth considered is 
sufficiently narrow, transmittance losses may be compensated 
for by using the average transmittance values for each while 
Still producing accurate results. The narrow bandwidth was 
the premise for the development of transmittance modelling 
in Chapter IV. 

Sensitivity of an image forming scanner is expressed by 
its noise equivalent temperature (NET). This parameter 
refers to the temperature difference AT within a scene 
element (with a given emissivity and average temperature of | 
e.g. 00০1) which can produce a change in the electrical 
signal level equal to the rms system noise. The noise in a 


well designed system is mostly detector noise. By 
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incorporating the atmospheric transmittance into equation 
6.8 and solving for AT, the NET is obtained as: 

NET - 4F VAf in °K (Eq. 6.9) 

৬0০0৮6০7548 

This an expression of the sensitivity of the detector system 
as a whole. It makes the assumption that the scene area 
over which the temperature difference AT occurs is much 
larger than the field of view. 

From the design aspect a major emphasis must be placed 
on the value of D* for a given detector since its value may 
vary as much as four orders of magnitude among common 
detector materials [Ref 2: Table 7.1]. In addition, the D* 
of a particular detector will vary considerably from the 
published data owing to such factors as manufacturing 
process, electrical bias, cooling, loading, aging and 
aerodynamic heating of the optical window. 

The electrical bandwidth (Af) of the system is 
determined by the scan rate. It is proportional to the 
reciprocal of the dwell time of an image point on the 
scanning detector surface. The amplifiers connected to the 
detector can be optimized to obtain a maximum signal-to- 
noise ratio, especially if the characteristics of the signal 
and noise are known. In this fashion a number of 
target/object discrimination techniques may be devised. 


Optimum filter design is extremely difficult to perform 
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using only a few components, though a simple band-pass 

filter with a flat response and a low and high frequency 

cutoff is often adequate. The very high signal to noise 

ratios realized by using matched filter (correlation?) 

techniques for processing are virtually impossible due to 

the passive nature of the systems, but by tweaking the high 

and low frequency roll-off a much improved background 
control and target definition may be accomplished with only 
a slight reduction of the S/N ratio. 

If a total field of view of 2 steradians is scanned at a 
constant rate without redundancy and with an instantaneous 
field of view W steradians, then the scan is completed in a 
frame time of ts. The dwell time t4 of a point on the 
detector is then ts-wW/0. As the total field of view is 
increased (either in azimuth, elevation or both) linearity 
in scan rate is difficult to achieve [Ref 4:pp 7351. The 
off-center loss in dwell time must be compensated for by a 
scan efficiency factor fg. As the scan efficiency 
decreases, an increased electrical bandwidth is required. 
In this case the optimum bandwidth required is 
approximately: 


Afra = SMEs in Hz (Eq. 6.10) 
4W ta 


Actually this refers to the upper cutoff frequency since it 


1s not practical to design for dc response. The lower 
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cutoff frequency can usually be designed to handle down to a. 
few Hertz which allows reproduction of at least most of the 
low frequency features in the scene. If a linear (vertical 
usually) array of detectors is used, with N elements in the 
array, then Afr, can be reduced proportionately. 

Equivalent noise bandwidth Af will usually be larger 
than required bandwidth Afr because of the inherent 1/f 
noise, amplifier noise, and R-C filter roll-off. In terms 
of scan rate, the noise equivalent temperature of a scan 
System is: 


πη... 
D.D“W, N 


NET = 2 


411 52 
ts 


W 6 6ے‎ 




















| 2 Sc: | in °K (Eq. 6.11) 


By breaking equation 6.11 into the above groups the 
corresponding components may be individually analyzed. The 
first group represents system performance requirements. The 
second group represents mostly the physical components 
available at design inception, and the third group 
represents efficiency parameters of materials used and the 
atmosphere. 

The principles of detectors and line scanning systems 
are covered extensively in References [23 and L31. The 
basics presented thus far will serve as a conceptual 
approach to modelling seeker systems, and will serve as a 
baseline for devising discrimination techniques in the 


modelling process. Based on discussion thus far no 
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committment has been made with regards to any particular 
system, but either conceived or existing seekers which use 
the line scanner techniques operate under the previously 


explained principles. 


fue Filtering for Discrimination 


There are three filtering methods which may be used 
in seeker system designs. These are: 1) Spectral, 2) 
Spatial, and 3) Temporal. Each method has its physical 
limitations for implementation in discrimination as will be 
discussed shortly. In contrast, the incorporation of 
discrimination into IR seekers leads to the necessity to 
analyze present countermeasures for possible enhancement to 


subjugate any of these discrimination techniques. 


a. Spectral Filtering 


Since background radiant intensity is a primary 
consideration to seeker design, methods to reduce it are a 
necessity. Background radiant intensity typically has peaks 
(maxima) at approximately 0.5 and 10 micrometers caused by 
solar reflection and thermal self-emission respectively. 

Around 3 microns a region of minimum background radiance ۱ 
exists. Coincidentally, between the two maxima, most hot 
targets produce their peak radiant intensity. The simplest 


method of detection is to examine typical target spectral 
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characteristics, select a detector most appropriate for 
operation at the target’s peak spectral points, and then use 
filtering to prevent the background radiant energy from 
entering the system. If the target of interest exhibits 
significant wavelength-dependent emissions, one method of 
discrimination could be the use of notch (2-color spectral 
bandpass) filters to select specific wavebancds. If the 
radiant energy in the passbands produced by the decoy is 
then less than that produced by the target (W/sr), the decoy 
nay be rendered ineffective in the seduction role. Thus 
much of the design of a decoy revolves around matching its 
spectral radiance to that of the ship it is being designed 
to protect. 

Most IR decoys are omni-directional by nature, 
and exhibit the same spectral signature in all directions 
(disregarding wind or sea-surftace effects if floating). 
Target signatures are aspect dependent (intensity) as well 
as spectrally dependent. As the aspect changes, so does the 
intensity of certain spectral components. The inherent 
variability of target signatures in this respect inakes it 
exceedingly difficult to establish discrimination techniques 
since look-up tables and matched filter methods cannot be 


used. 
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b. Temporal Filtering 


Spatial frequencies produce time varying signals 
with a frequency content equal to the spatial frequency 
times the angular scan rate. In processing then, low 
frequency rejection and high frequency peaking can be used 
to eliminate background clutter. Only if the target and 
background signals can be adequately described can a spatial 
filter maximize the ratio of target signal to background 
Signal. If multi-element detectors are used, simpler filter 
techniques are possible. The instantaneous field would be 
made as small as practical, and the upper frequency cutoff 
of the electrical filter is selected to correspond to the 
dwell time of a point on the detector. The low frequency 
cutoff is usually set at 1/3 t5 1/4 of the upper cutoff 
frequency to Block the Background radiance without allowing 
saturation at t^ other end of the bandpass. This filtering 
technique 1s most often used to decrease the probability of 
false alarms through scan-to-scan comparisons and has little 


merit in discerning targets from decoys. 


pM Spatial Filtering 


A more promising approach to discrimination by 
the seeker may be to process spatial characteristics of the 
scene. When the range from the seeker to the target is 


large, the objects in the scene will appear as point 
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sources. As this range decreases, however, large objects 
form a distribution within the scene (both horizontal and 
vertical) while small objects retain their point source 
characteristics. By keeping the instantaneous field of view 
small, integral (in spatial sense) processing could be used 
to discern large distributed objects from small ones. It 
then becomes the task of the decoy designer to not only make 
the decoy more intense spectrally than the ship, but to make 
the signature more widely distributed to simulate large 
targets. If the decoy is suspended in the air, for example 
by parachute (and somewhat less so if the decoy is 
waterborne), the wind will tend to distribute the heat over 
a wider area; but variable wind conditions and short hang 
times of parachute-deployed decoys make this method of 
deployment less attractive than water-borne placement. Some 
of the effects of wind for signature expansion are still 
present with water-borne decoys. In addition to the torch 
emission above the sea surface there exists a reflective 
component off the sea surface referred to as the "Halo 
effect" (Ref. 4]. This effect was originally noticed by 
personnel from Naval Research Laboratory taking infrared 
measurements of the USS Ticonderoga’s signature in both the 
5-95 and 8-14 micron bands. The halo effect may readily be 
seen in Figure 6-1. Upon investigation it was established 
that the sea foam caused by a white-capped sea surface has 


very high reflectance at low grazing angles in both infrared 


79 


Figure 6-1. 
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Haloing caused by sea-surface 
reflection around USS Ticonderoga. 
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bands. The white-caps may be caused by wind, ship movement, - 
and other surface turbulence. When viewed from certain 
angles as in Figure 6-2, energy from the source which was 
not originally emitted in the direction of the seeker 
(observer) is reflected from the sea surface toward the 
seeker, adding to the original signature radiant density. 

From the decoy design aspect this effect may be 
capitalized upon by incorporating a means of disturbing 
(foaming) the sea surface over a wide area around a water- 
borne decoy beneath its burning source. By floating the 
decoy at the proper height above the foamed sea surface a 
very large reflective signature may be achieved in addition 
to the intense point source of the torch itself. 

Foaming may be accomplished in many ways, but 
primary interest must be in distributing over as large an 
area as possible. Since it is difficult to foam the sea 
surface over a wide area using a simple squirrel-cage fan or 
compressed air (too many cylinders) another prospect must be 
taken which considers weight and space restrictions of the 
decoy. Foaming of sea water is purely transient and very 
short in duration, so the above methods would probably be 
discarded upon testing. The alternate means suggested here 
is to use a small pneumatically driven suction pump with a 
mixing chamber to mix the incoming sea water with, for 
example, light-water or aqueous film-forming foam (AFFF). 


The mixture may then be sprayed from the decoy by fan or 
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nozzles over a wide area of the sea surface to create a 
highly reflective layer (using sea surface turbulence for 
wide angular reflection) and correspondingly large 
signature. The components in this case could be very 
inexpensive. Additionally, this technique would work in all 
infrared bands and, because of the difficulty in devising 
passive-only discrimination methods, may be the optimum 
offboard countermeasure to IK seekers. The methods 
discussed here have not been attempted, but these and many 
other approaches will need to be addressed in order to meet 


the IR threat with inexpensive offboard countermeasures. 
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Wind=5 


Wind-10 


Wind-15 


Wind-20 


TRANSMITTANCE IN 


ΧΧΧ 


TABLE 


A. i. ay 


3.5 UM BAND - TROPICAL PROFILE 
LOWTRAN DATA 


ΧΧΧ 


Visibility in Kilometers 


30 


20 


15 


. 16% 


. 50% 
. 30% 
. 10% 
. 30% 
«90. 


. 76% 
. 84% 


21% 


. 41% 
. 83% 
. 90% 
. 80% 
‘67% 


. 90% 
‘61% 
. 81% 
. 85% 
. 21% 
. 38% 
. 43% 
. 38% 


. 87X 
. 50% 
. 62% 
. 59X 
. 92% 
. 
- 27% 
. 


.75X 


3.70% 


. 47% 
. 83% 
. 88% 
. 73% 
. 88% 
10.84% 
4.21% 
2.67% 


- 39% 
.51X 
.37X 
. 03% 
. 19% 
10.30% 
3. 90% 
2.45% 


. 35% 
. 36% 
. 12% 
. 71% 
19.87% 
10.03% 
3.77% 
2.33% 
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93% 


. 44% 


. 49% 
. 92% 
. 70% 
. 68% 
. 60% 
. 70% 
.61% 
. 87% 


.37X 
. 46% 
. 98% 
. 80% 
ZU 
.27X 
. 43% 
. 76% 


. 32% 
. 24% 
. 65% 
. 40% 
. 59% 
. 08% 
- 35% 
. 71% 


.11% 


.55% 


. 52% 
. 19% 
. 10% 
. 13X 
. 83% 
. 98% 
. 61x 
. 28% 


. 36% 
. 59% 
228 
. 17% 
. 22% 
. 68% 
.52X 
.23X 


. 28% 
. 31% 
. 82% 
. 74% 
. 95% 
. 55% 
. 48% 
.21X 


.61X 


. 37X 


. 49% 
. 68% 
. 33% 
. 69% 
. 07% 
. 026 


. 26% 
.11% 
. 00% 
. 25% 
. 17% 
. 64% 
. 06% 
. 02% 


— TABLE Α.].Ὀὺ - 
TRANSMITTANCE IN 3.5 UM BAND - TROPICAL PROFILE 
KKK EMPIRICAL “xk 
Range !: Visibility in Kilometers 
in Km | 40 30 20 15 10 

Wind= 0 0.2 : 86.63% 86.58% 86.57% 86.49% 86.33% 
Wind= 0 1.0: 81.37% 81.13% 81.08% 80.72% 79.94% 
Wind= 0 2.0: 75.25% 74.80% 74.70% 74.05% 72.63% 
Wind= 0 4.0 : 64.34% 63.59% 63.41% 62.31% 59.94% 
Wind= 0 8.0 : 47.04% 45.95% 45.69% 44.12% 40.82% 
Wind= 012.0 : 34.40% 33.20% 32.93% 31.24% 27.80% 
Wind= 0 18.0: 21.50% 20.39% 20.14% 18.61% 15.63% 
Wind= 0 21.0 : 17.00% 15.98% 15.75% 14.36% 11.72% 
Wind= 5 0.2 : 85.75% 85.14% 84.16% 83.06% 80.76% 
Wind= 5 1.0: 77.30% 74.60% 70.42% 65.94% 57.28% 
Wind- 5 2.0 : 67.91% 63.23% 56.35% 49.41% 37.29% 
Wind= 5 4.0 : 52.40% 45.44% 36.08% 27.74% 15.80% 
Wind- 5 8.0 : 31.21% 23.46% 14.79% 8.74% 2.84% 
Wind- 5 12.0 : 18.58% 12.11% 6.07% 2.76% 0.51% 
Wind- 5 18.0 ¦ 8.54% 4.49% 1.59% 0.49% 0.04% 
Wind= 5 21.0 : 5.79% 2.74% 0.82% 0.21% 0.01% 
Wind=-10 0.2 : 85.55% 84.88% 83.82% 82.63% 80.13% 
Wind=10 1.0: 76.39% 73.48% 69.00% 64.22% 55.08% 
Wind-10 2.0 : 66.32% 61.35% 54.10% 46.87% 34.47% 
Wind=10 4.0 : 49.98% 42.77% 33.26% 24.96% 13.50% 
Wind-10 8.0 : 28.38% 20.79% 12.57% 7.08% 2.07% 
Wind=10 12.0 : 16.12% 10.10% 4.75% 2.01% 0.32% 
Wind=10 18.0 : 6.90% 3.42% 1.10% 0.30% 0.02% 
Wind-10 21.0 ; 4.51% 1.99% 0.53% 0.12% 0.00% 
Wind=15 0.2 : 85.50% 84.82% 83.74% 82.53% 79.99% 
Wind-15 1.0: 76.19% 73.23% 68.68% 63.84% 54.60% 
Wind-15 2.0 : 65.97% 60.94% 53.61% 46.32% 33.87% 
Wind-15 4.0 : 49.45% 42.19% 32.66% 24.38% 13.04% 
Wind-15 8.0 : 27.79% 20.23% 12.12% 6.75% 1.93% 
Wind=15 12.0 : 15.61% 9.70% 4.50% 1.87% 0.29% 
Wind-15 18.0 : 6.58% 3.22% 1.02% 0.27% 0.02% 
Wind-15 21.0 ; 4.27% 1.86% 0.48% 0.10% 0.00% 
Wind=20 0.2 : 85.49% 84.81% 83.73% 82.51% 79.96% 
Wind=20 1.0 : 76.15% 73.17% 68.61% 63.76% 54.49% 
Wind=20 2.0 : 65.89% 60.84% 53.50% 46.19% 33.74% 
Wind=20 4.0 : 49.33% 42.07% 32.52% 24.25% 12.94% 
Wind-20 8.0 : 27.66% 20.11% 12.02% 6.68% 1.90% 
Wind=20 12.0 : 15.50% 9.61% 4.44% 1.84% 0.28% 
Wind-20 18.0 : 6.51% 3.18% 1.00% 0.27% 0.02% 
Wind=20 21.0 | 4.22% 1.83% 0.47% 0.10% 0.00% 
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- TABLE A.2 - 


TRANSMITTANCE IN 3.5 UM BAND - MIDLATITUDE SUMMER PROFILE 


184 LOWTRAN DATA 88184 
Range ! Visibility in Kilometers 
in Km : 40 30 20 15 10 5 
Wind=0 0.20 : 94.16% 94.12% 94.06% 94.00% 93.87% 93.49% 
1.00 | 84.47% 84.33% 84.05% 83.76% 83.20% 81.53% 
2.00 | 76.79% 76.53% 76.02% 75.51% 74.49% 71.53% 
4.00 : 65.59% 65.16% 64.28% 63.42% 61.72% 56.92% 
8.00 : 50.26% 49.59% 48.27% 46.98% 44.50% 37.84% 
12.00 ! 39.82% 39.03% 37.48% 35.99% 33.17% 26.02% 
18.00 : 28.86% 28.01% 26.36% 24.79% 21.94% 15.24% 
21.00 : 24.90% 24.04% 22.40% 20.86% 18.09% 11.82% 
Wind=5 0.20 : 93.14% 92.72% 91.89% 91.06% 89.43% 84.72% 
1.00 : 80.01% 78.24% 74.80% 71.48% 65.29% 49.83% 
2.00 : 68.89% 65.88% 60.21% 54.99% 45.88% 26.72% 
4.00 |: 52.78% 48.28% 40.33% 33.64% 23.41% 7.94% 
8.00 ؛‎ 32.55% 27.23% 19.00% 13.22% 6.40% 0.74% 
12.00 : 20.75% 15.88% 9.26% 5.37% 1.81% 0.07% 
18.00 : 10.86% 7.27% 3.23% 1.43% 0.28% 0.00% 
21.00 : 7.96% 4.98% 1.94% 0.75% 0.11% 0.00% 
Wind-10 0.20 | 92.93% 92.44% 91.45% 90.47% 88.53% 83.01% 
1.00 : 79.11% 77.04% 73.02% 69.18% 62.10% 45.00% 
2.00 : 67.35% 63.86% 57.38% 51.50% 41.50% 21.79% 
4.00 : 50.46% 45.37% 36.62% 29.50% 19.15% 5.28% 
8.00 : 29.74% 24.05% 15.66% 10.16% 4.29% 0.33% 
12.00 : 18.13% 13.18% 6.93% 3.62% 0.99% 0.02% 
18.00 ¦ 8.87% 5.49% 2.09% 0.79% 0.11% 0.00% 
21.00 : 6.28% 3.60% 1.17% 0.38% 0.04% 0.00% 
Wind-15 0.20 : 92.87% 92.36% 91.33% 90.30% 88.29% 82.55% 
1.00 : 78.87% 76.71% 72.54% 68.55% 61.25% 43.76% 
2.00 : 66.94% 63.32% 56.62% 50.58% 40.37% 20.60% 
4.00 : 49.84% 44.60% 35.66% 28.45% 18.13% 4.72% 
8.00 : 29.02% 23.24% 14.86% 9.46% 3.84% 0.26% 
12.00 : 17.47% 12.52% 6.40% 3.25% 0.84% 0.01% 
18.00 : 8.39% 5.09% 1.86% 0.67% 0.09% 0.00% 
21.00 ¦ 5.89% 3.29% 1.02% 0.31% 0.03% 0.00% 
Wind-20 0.20 : 92.84% 92.32% 91.27% 90.23% 88.18% 82.33% 
1.00 : 78.75% 76.56% 72.31% 68.27% 60.85% 43.19% 
2.00 : 66.74% 63.07% 56.27% 50.15% 39.85% 20.07% 
4.00 : 49.55% 44.25% 35.22% 27.98% 17.66% 4.48% 
8.00 |: 28.68% 22.87% 14.49% 9.14% 3.65% 0.23% 
12.00 ؛‎ 17.۱۶7۹ ae. 22% 6.16% 3.09% 0.78% 0.01% 
18.00 ; 8.17% 4.91% 1.76% 0.62% 0.08% 0.00% 
21900 0 ٤ 3.15% 0.95% 0.28% 0.03% 0.00% 
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— TABLE A.2.b - 


TRANSMITTANCE IN 3.5 UM BAND —- MIDLATITUDE SUMMER PROFILE 
84858 EMPIRICAL xxx 


Range : Visibility in Kilometers 
in Km : 40 30 20 15 10 5 
Wind= 0 0.2 : 86.48% 86.39% 86.35% 86.26% 86.09% 85.79% 
Wind= 0 1.0: 82.52% 82.07% 81.88% 81.47% 80.65% 79.28% 
Wind= 0 2.0 ! 77.82% 76.99% 76.62% 75.86% 74.34% 71.83% 
Wind= 0 4.0 : 69.21% 67.73% 67.09% 65.77% 63.16% 58.96% 
Wind= 0 8.0 : 54.74% 52.43% 51.45% 49.44% 45.60% 39.73% 
Wind= 0 12.0 : 43.30% 40.59% 39.45% 37.16% 32.92% 26.77% 
Wind= 0 18.0 : 30.46% 27.64% 26.49% 24.22% 20.19% 14.81% 
Wind= 0 21.0 : 25.55% 22.81% 21.70% 19.55% 15.81% 11.01% 
Wind= 5 0.2 : 85.70% 85.07% 84.11% 83.02% 80.69% 74.23% 
Wind= 5 1.0: 78.85% 76.01% 71.83% 67.29% 58.34% 38.46% 
Wind= 5 2.0 : 71.06% 66.03% 58.97% 51.74% 38.89x 16.90% 
Wind= 5 4.0 : 57.70% 49.83% 39.74% 30.60% 17.29% 3.26% 
Wind- 5 8.0 : 38.05% 28.38% 18.05% 10.70% 3.42% 0.12% 
Wind- 5 12.0 : 25.10% 16.16% 8.20% 3.74% 0.67% 0.00% 
Wind= 5 18.0 : 13.44% 6.95% 2.51% 0.77% 0.06% 0.00% 
Wind- 5 21.0 : 9.84% 4.55X 1.39% 0.35% 0.02% 0.00% 
Wind-10 0.2 : 85.49% 84.79% 83.73% 82.52% 79.93% 72.83% 
Wind=10 1.0: 77.90% 74.78% 70.20% 65.26% 55.65% 34.95% 
Wind-10 2.0 : 69.36% 63.90% 56.32% 48.67% 35.39% 13.96% 
Wind-10 4.0 : 54.97% 46.67% 36.25% 27.08% 14.32% 2.23% 
Wind=10 8.0 ! 34.54% 24.90% 15.02% 8.38% 2.34% 0.06% 
Wind=10 12.0 : 21.70% 13.28% 6.22% 2.59% 0.38% 0.00% 
Wind=10 18.0 : 10.81% 5.17% 1.66% 0.45% 0.03% 0.00% 
Wind-10 21.0 : 7.63% 3.23% 0.86% 0.19% 0.01% 0.00% 
Wind=15 0.2 : 85.44% 84.72% 83.63% 82.39% 79.73% 72.47% 
Wind=15 1.0: 77.66% 74.46% 69.78% 64.75% 54.97% 34.09% 
Wind-15 2.0 1: 68.92% 63.36% 55.65% 47.91% 34.54% 13.28% 
Wind-15 4.0 : 54.29% 45.89% 35.40% 26.23% 13.63% 2.02% 
Wind=15 8.0 | 33.68% 24.06% 14.32% 7.86% 2.12% 0.05% 
Wind=15 12.0 : 20.90% 12.62% 5.79% 2.36% 0.33% 0.00% 
Wind-15 18.0 : 10.21% 4.79% 1.49% 0.39% 0.02% 0.00% 
Wind=15 21.0 : 7.14% 2.95% 0.76% 0.16% 0.01% 0.00% 
Wind-20 0.2 ؛‎ 85.42% 84.70% 83.60% 82.35% 79.68% 72.37% 
Wind=20 1.0: 77.59% 74.38% 69.68% 64.61% 54.80% 33.87% 
Wind=20 2.0 ! 68.81% 63.23% 55.48% 47.71% 34.32% 13.11% 
Wind=20 4.0 : 54.11% 45.68% 35.18% 26.02% 13.46% 1. 96% 
Wind=20 8.0 : 33.46% 23.85% 14.15% 7.74% 2.07% 0.04% 
Wind=20 12.0 : 20.70% 12.45% 5.69% 2.30% 0.32% 0.00% 
Wind=20 18.0 : 10.06% 4.70% 1.45% 0.37% 0.02% 0.00% 
Wind=20 21.0 ; 7.02% 2.89% 0.73% 0.15% 0.00% 0.00% 
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— TABLE SRA, 3. a | 


TRANSMITTANCE IN 3.5 UM BAND —- MIDLATITUDE WINTER PROFILE 
ΧΧΧ LOWTRAN DATA ΧΧΧ 


Range Visibility in Kilometers 
in Km 40 30 20 15 10 5 
Wind=0 0.2 ¦ 97.55% 97.52% 97.45% 97.39% 97.26% 96.86% 
1.00 : 92.46% 92.30% 91.99% 91.68% 91.06% 89.22% 
2.00 : 88.34% 88.04% 87.45% 86.85% 85.68% 82.26% 
4.00 : 82.24% 81.69% 80.59% 79.50% 77.36% 71.32% 
8.00 : 73.25% 72.28% 70.35% 68.45% 64.82% 55.09% 
12.00 . 66.20% 64.88% 62.30% 59.80% 55.10% 43.18% 
18.00 : 57.51% 55.80% 52.51% 49.38% 43.68% 30.29% 
21.00 : 53.76% 51.90% 48.34% 45.00% 39.00% 25.45% 
Wind=5 0.2 : 96.50% 96.07% 95.21% 94.35% 92.66% 87.78% 
1.00 ¦ 87.57% 85.64% 81.88% 78.25% 71.47% 54.55% 
2.00 |: 79.25% 75.79% 69.27% 63.27% 52.78% 30.75% 
4.00 : 66.19% 60.54% 50.58% 42.18% 29.36% 9.96% 
8.00 : 47.45% 39.69% 27.70% 19.27% 9.34% 1.08% 
12.00 : 34.51% 26.40% 15.40% 8.93% 3.01% 0.12% 
18.00 : 21.65% 14.49% 6.45% 2.85% 0.56% 0.00% 
21.00 : 17.19% 10.76% 4.19% 1.62% 0.24% 0.00% 


Wind=10 0.2 96.28% 95.77% 94.75% 93.73% 91.72% 85.99% 


1.00 : 86.58% 84.31% 79.90% 75.69% 67.94% 49.21% 
2.00 . 77.47% 73.45% 65.98% 59.21% 47.69% 25.02% 
4.00 : 63.24% 56.86% 45.87% 36.94% 23.97% 6.60% 
8.00 | 43.32% 35.01% 22.79% 14.78% 6.22% 0.47% 
12.00 : 30.11% 21.87% 11.49% 6.00% 1.64% 0.03% 
18.00 : 17.64% 10.92% 4.16% 1.57% 0.22% 0.00% 
21.00 : 13.54% 7.74% 255 X 0.81% 0.08% 0.00% 


Wind=15 0. 2 96.22% 95.69% 94.62% 93.56% 91.47% 85.51% 


47.83% 67.00% 75.01% 79.37% 83.95% 86.31% : 1.00 
23.64% 46.38% 58.13% 65.10% 72.82% 76.98% | 2.00 
5.89% 22.67% 35.62% 44.66% 55.88% 62.46% | 4.00 
0.38% 5.57% 13.74% 21.60% 33.82% 42.25% ؛ 8.00 
0.02% 1.39% 5.38% 10.60% 20.77% 29.00% : 12.00 
0.00% 0.17% 1.33% 3.69% 10.11% 16.68% ؛ 18.00 
0.00% 0.06% 0.66% 2.18% 7.07% 12.68% ؛ 21.00 


Wind=20 0.2 96.19% 95.65% 94.56% 93.48% 91.35% 85.28% 


1.00 : 86.19% 83.78% 79.12% 74.69% 66.56% 47.20% 
2.00 , 76.76% 72.52% 64.69% 57.64% 45.78% 23.02% 
4.00 : 62.09% 55.43% 44.10% 35.01% 22.09% 5.59% 
8.00 : 41.76% 33.28% 21.07% 13.28% 5.28% 0.34% 
12.00 : 28.49% 20.27% 10.21% 5.11% 1.28% 0.02% 
18.00 : 16.24% 9.75% 3.48% 1.23% 0.16% 0.00% 
21.00 : 12.30% 6.78% 2.04% 0.61% 0.05% 0.00% 


88 


- TABLE 4.3.b - 


TRANSMITTANCE IN 3.5 UM BAND _ MIDLATITUDE WINTER PROFILE 
184 EMPIRICAL 184 


Range : Visibility in Kilometers 
in Km : 40 30 20 15 10 5 
Wind= 0 0.2 : 93.30% 93.24% 93.23% 93.15% 92.97% 92.65% 
Wind= 0 1.0: 91.21% 90.93% 90.87% 90.47% 89.59% 88.08% 
Wind= 0 2.0 : 88.66% 88.12% 87.99% 87.22% 85.54% 82.68% 
Wind= 0 4.0 : 83.77% 82.75% 82.52% 81.08% 77.99% 72.85% 
Wind- O 8.0 : 74.78% 72.98% 72.57% 70.06% 64.82% 56.56% 
Wind= 0 12.0 : 66.76% 64.37% 63.83% 60.54% 53.87% 43.92% 
Wind= 0 18.0 : 56.32% 53.31% 52.64% 48.63% 40.82% 30.05% 
Wind= 0 21.0 : 51.72% 48.52% 47.81% 43.59% 35.53% 24.85% 
Wind= 5 0.2 : 92.32% 91.65% 90.61% 89.09% 85.29% 73.25% 
Wind- 5 1.0 : 86.51% 83.44% 78.80% 72.43% 58.24% 27.21% 
Wind= 5 2.0 : 79.75% 74.19% 66.18% 55.90% 36.15% 7.89% 
Wind= 5 4.0 : 67.79% 58.67% 46.67% 33.31% 13.93% 0.66% 
Wind= 5 8.0 : 48.97% 36.68% 23.22% 11.82% 2.07% 0.00% 
Wind= 5 12.0 : 35.38% 22.94% 11.55% 4.20% 0.31% 0.00% 
Wind= 5 18.0 : 21.73% 11.34% 4.05% 0.89% 0.02% 0.00% 
Wind= 5 21.0 : 17.03% 7.97% 2.40% 0.41% 0.00% 0.00% 
Wind-10 0.2 : 92.09% 91.33% 90.14% 88.41% 84.08% 70.58% 
Wind-10 1.0: 85.46% 81.98% 76.77% 69.67% 54.22% 22.60% 
Wind-10 2.0 : 77.83% 71.63% 62.81% 51.73% 31.34% 5.44% 
Wind=10 4.0 : 64.57% 54.68% 42.04% 28.52% 10.46% 0.32% 
Wind-10 8.0 : 44.43% 31.87% 18.83% 8.67% 1.17% 0.00% 
Wind=10 12.0 : 30.57% 18.57% 8.44% 2.63% 0.13% 0.00% 
Wind-10 18.0 : 17.45% 8.26% 2.53% 0.44% 0.00% 0.00% 
Wind=-10 21.0: 13.18% 5.51% 1.39% 0.18% 0.00% 0.00% 
Wind-15 0.2 : 92.04% 91.26% 90.03% 88.25% 83.81% 69.99% 
Wind-15 1.0: 85.22% 81.65% 76.31% 69.05% 53.34% 21.66% 
Wind-15 2.0 : 77.40% 71.06% 62.06% 50.82% 30.33% 5.00% 
Wind-15 4.0 : 63.85% 53.81% 41.04% 27.53% 9.80% 0.27% 
Wind-15 8.0 : 43.45% 30.86% 17.95% 8.07% 1.02% 0.00% 
Wind=15 12.0 : 29.57% 17.70% 7.85% 2.37% 0.11% 0.00% 
Wind=15 18.0 : 16.60% 7.69% 2.27% 0.38% 0.00% 0.00% 
Wind=15 21.0 : 12.43% 5.06% 1.22% 0.15% 0.00% 0.00% 
Wind=20 0.2 : 92.03% 91.24% 90.01% 88.21% 83.75% 69.85. 
Wind=20 1.0: 85.17% 81.58% 76.20% 68.91% 53.15% 21.45% 
Wind=20 2.0 : 77.30% 70.92% 61.89% 50.61% 30.10% 4.91% 
Wind=20 4.0 : 63.69% 53.61% 40.82% 27.30% 9.66% 0.26% 
Wind-20 8.0 : 43.23% 30.63% 17.76% 7.94% 0.99% 0.00% 
Wind=20 12.0 : 29.34% 17.50% 7.73% 2.31% 0.10% 0.00% 
Wind=20 18.0 : 16.41% 7.56% 2.22% 0.36% 0.00% 0.00% 
Wind=20 21.0 : 12.27% 4.97% 1.19% 0.14% 0.00% 0.00% 
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Wind-0 


Wind-5 


Wind-10 


Wind-15 


Wind-20 


— TABLE A.4.a - 


TRANSMITTANCE IN 3.5 UM BAND - SUB-ARCTIC SUMMER PROFILE 
114 LOWTRAN DATA 11. 


Range ا‎ Visibility in Kilometers 
in Km : 40 Km 30 Km 20 Km 15 Km 10 Km 5 Km 
0.20 : 95.52% 95.49% 95.43% 95.36% 95.23% 94.85% 
1.00 ; 87.64% 87.49% 87.20% 86.91% 86.32% 84.60% 
2.00 : 81.41% 81.14% 80.60% 80.05% 78.98% 75.86% 
4.00 : 72.23% 71.75% 70.79% 69.84% 67.98% 62.72% 
8.00 : 59.00% 58.22% 56.68% 55.17% 52.26% 44.48% 
12.00 : 49.49% 48.51% 46.60% 44.74% 41.26% 32.40% 
18.00 ; 38.90% 37.75% 35.54% 33.44% 29.61% 20.60% 
21.00 ؛‎ 34.70% 33.50% 31.23% 29.09% 25.24% 16.53% 
0.20 . 94.49% 94.07% 93.23% 92.38% 90.73% 85.95% 
1.00 : 83.00% 81.17% 77.60% 74.16% 67.73% 51.69% 
2.00 : 73.02% 69.84% 63.83% 58.29% 48.63% 28.32% 
4.00 ؛‎ 58.11% 53.15% 44.40% 37.03% 25.77% 8.74% 
8.00 : 38.19% 31.95% 22.30% ΙΟ 5 7.51% 0.86% 
12.00 ؛‎ 25.78% 19.72% 11.50% 6.67% 2.25% 0.09% 
18.00 ؛‎ 14.62% 9.79% 4.36% 1.92% 0.38% 0.00% 
21.00 ; 11.08% 6.93% 2.70% 1.04% 0.16% 0.00% 
0.20 | 94.28% 93.78% 92.78% 91.79% 89.83% 84.23% 
1.00 : 82.08% 79.93% 75.77% 71.79% 64.45% 46.73% 
2.00 : 71.42% 67.72% 60.85% 54.63% 44.03% 23.15% 
4.00 : 55.58% 49.98% 40.36% 32.52% 21.13% 5.84% 
8.00 | 34.94% 28.26% 18.42% 11.96% 5.05% 0.39% 
12.00 ; 22.55% 16.40% 8.63% 4.52% 1.24% 0.03% 
18.00 : 11.97% 7.42% 2.83% 1.07% 0.15% 0.00% 
21.00 , 8.77% 5.02% 1.63% 0.53% 0.05% 0.00% 
0.20 : 94.22% 93.70% 92.66% 91.62% 89.58% 83.77% 
1.00 : 81.83% 79.60% 75.28% 71.15% 63.58% 45.46% 
2.00 ؛‎ 70.98% 67.16% 60.06% 53.66% 42.85% 21.90% 
4.00 : 54.91% 49.15% 39.31% 31.38% 20.01% 5.23% 
8.00 : 34.10% 27.32% 17.48% ΠΣ 4.53% 0.31% 
12.00 : 21.74% 15.59% 7.98% 4.06% 1.05% 0.02% 
18.00 ; 11.33% 6.88% 252% 0.91% 0.12% 0.00% 
21.00 ; 8.23X 4.60% 1.42% 0.44% 0.04% 0.00% 
0.20 ؛‎ 94.20% 93.67% 92.60% 91.55% 89.47% 83.55% 
1.00 : 81.72% 79.44% 75.05% 70.86% 63.18% 44.87% 
2.00 : 70.78% 66.89% 59.70% 53.22% 42.31% 21.34% 
4.00 : 54.60% 48.76% 38.84% 30.86% 19.51% 4.96% 
8.00 : 33.71% 26.89% 17.06% 10.77% 4.30% 0.28% 
12.00 : 20% 37% 15.253 7.69% 3.86% 0.97% 0.02% 
18.00 ; 11.04% 6.64% 2.38% 0.85% 0.11% 0.00% 
21.00 ; 7.98% 4.41% 1.34% 0.40% 0.04% 0.00% 
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TRANSMITTANCE IN 3.5 UM BAND —- SUB-ARCTIC SUMMER 
EMPIRICAL 
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Visibility in Kilometers 
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TRANSMITTANCE IN 3.5 UM BAND - SUB-ARCTIC 
LOWTRAN DATA 


ΧΧΧ 


— TABLE 
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Visibility in Kilometers 
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— TABLE A.5.b - 
TRANSMITTANCE IN 3.5 UM BAND - SUBARCTIC WINTER PROFILE 
KKK EMPIRICAL 1184 

Range | Visibility in Kilometers 

in Km . 40 30 20 15 10 5 
Wind= 0 0.2 : 96.02% 95.96% 95.95% 95.86% 95.68% 95.34% 
Wind= 0 1.0: 94.69% 94.40% 94.33% 93.91% 92.98% 91.39% 
Wind= 0 2.0 : 93.05% 92.47% 92.34% 91.52% 89.73% 86.67% 
Wind= 0 4.0 | 89.85% 88.74% 88.49% 86.92% 83.55% 77.96% 
Wind= 0 8.0: 83.78% 81.73% 81.26% 78.41% 72.44% 63.08% 
Wind= 0 12.0 : 78.12% 75.27% 74.63% 70.73% 62.81% 51.03% 
Wind= 0 18.0 : 70.34% 66.52% 65.67% 60.59% 50.71% 37.14% 
Wind= 0 21.0 : 66.74% 62.54% 61.61% 56.09% 45.57% 31.68% 
Wind- 5 0.2 : 95.02% 94.28% 93.17% 91.28% 86.18% 69.38% 
Wind= 5 1.0 : 89.86% 86.39% 81.42% 73.51% 55.13% 18.65% 
Wind- 5 2.0 : 83.80% 77.46% 68.80% 56.08% 31.54% 3.61% 
Wind= 5 4.0 | 72.87% 62.26% 49.12% 32.64% 10.32% 0.14% 
Wind- 5 8.0 : 55.11% 40.23% 25.04% 11.06% 1.11% 0.00% 
Wind= 5 12.0 : 41.68% 26.00% 12.76% 3.75% 0.12% 0.00% 
Wind= 5 18.0 |: 27.41% 13.50% 4.65% 0.74% 0.00% 0.00% 
Wind= 5 21.0 1 22.23% 9.73% 2.80% 0.33% 0.00% 0.00% 
Wind-10 0.2 : 94.78% 93.90% 92.60% 90.38% 84.44% 65.34% 
Wind=-10 1.0: 88.72% 84.69% 78.95% 69.96% 49.78% 13.81% 
Wind=10 2.0 !' 81.68% 74.43% 64.69% 50.80% 25.71% 1.98% 
Wind=-10 4.0 . 69.24% 57.49% 43.43% 26.78% 6.86% 0.04% 
Wind-10 8.0 : 49.76% 34.30% 19.57% 7.44% 0.49% 0.00% 
Wind=-10 12.0 : 35.76% 20.46% 8.82% 2.07% 0.03% 0.00% 
Wind=-10 18.0 : 21.78% 9.43% 2.67% 0.30% 0.00% 0.00% 
Wind-10 21.0 : 17.00% 6.40% 1.47% 0.12% 0.00% 0.00% 
Wind=15 0.2 ؛‎ 94.72% 93.81% 92.46% 90.17% 84.02% 64.39% 
Wind-15 1.0: 88.44% 84.28% 78.36% 69.13% 48.56% 12.84% 
Wind-15 2.0 : 81.18% 73.71% 63.73% 49.59% 24.47% 1.71% 
Wind-15 4.0 : 68.39% 56.38% 42.14% 25.52% 6.21% 0.03% 
Wind-15 8.0 |: 48.54% 32.99% 18.43% 6.76% 0.40% 0.00% 
Wind-15 12.0 | 34.45% 19.30% 8.06% 1.79% 0.03% 0.00% 
Wind=15 18.0 : 20.60% 8.64% 2.33% 0.24% 0.00% 0.00% 
Wind=15 21.0 : 15.93% 5.78% 1.25% 0.09% 0.00% 0.00% 
Wind=20 0.2 : 94.71% 93.79% 92.42% 90.11% 83.91% 64.16% 
Wind=20 1.0: 88.38% 84.18% 78.22% 68.92% 48.26% 12.61% 
Wind=20 2.0 : 81.06% 73.53% 63.49% 49.30% 24.17% 1.65% 
Wind=20 4.0 |: 68.18% 56.11% 41.84% 25.22% 6.06% 0.03% 
Wind-20 8.0 : 48.24% 32.68% 18.17% 6.60% 0.38% 0.00% 
Wind-20 12.0 : 34.14% 19.03% 7.89% 1.73% 0.02% 0.00% 
Wind-20 18.0 : 20.32% 8.46% 2.26% 0.23% 0.00% 0.00% 
Wind=20 21.0 : 15.67% 5.64% 1.21% 0.08% 0.00% 0.00% 
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— TABLE A.6.b - 
TRANSMITTANCE IN 3.5 UM BAND - STANDARD ATMOSPHERE PROFILE 
ΧΧΧ EMPIRICAL 114 
Range | Visibility in Kilometers 
in Km : 40 30 20 15 10 5 

Wind= 0 0.2 : 91.41% 91.36% 91.35% 91.28% 91.12% 90.84% 
Wind= 0 1.0: 88.78% 88.53% 88.48% 88.13% 87.37% 86.06% 
Wind= 0 2.0 : 85.59% 85.13% 85.02% 84.36% 82.91% 80.43% 
Wind= 0 4.0 : 79.56% 78.70% 78.50% 77.28% 74.65% 70.26% 
Wind= 0 8.0 : 68.73% 67.26% 66.92% 64.86% 60.52% 53.61% 
Wind= 0 12.0 : 59.38% 57.48% 57.06% 54.44% 49.06% 40.90% 
Wind= 0 18.0 : 47.69% 45.42% 44.91% 41.86% 35.81% 27.26% 
Wind= 0 21.0 : 42.74% 40.37% 39.85% 36.70% 30.60% 22.26% 
Wind- 5 0.2 : 90.45% 89.84% 88.86% 87.58% 84.55% 75.35% 
Wind= 5 1.0: 84.21% 81.41% 77.08% 71.67% 60.11% 33.78% 
Wind= 5 2.0 , 77.01% 71.98% 64.53% 55.79% 39.23% 12.39% 
Wind- 5 4.0 : 64.41% 56.26% 45.22% 33.80% 16.72% 1.67% 
Wind- 5 8.0 : 45.06% 34.38% 22.20% 12.40% 3.04% 0.03% 
Wind= 5 12.0 : 31.52% 21.01% 10.90% 4.55% 0.55% 0.00% 
Wind= 5 18.0 : 18.44% 10.03% 3.75% 1.01% 0.04% 0.00% 
Wind= 5 21.0 : 14.11% 6.93% 2.20% 0.48% 0.01% 0.00% 
Wind-10 0.2 ; 90.27% 89.60% 88.52% 87.10% 83.77% 73.72% 
Wind-10 1.0: 83.39% 80.33% 75.60% 69.74% 57.37% 30.28% 
Wind-10 2.0 : 75.53% 70.07% 62.07% 52.82% 35.75% 9.96% 
Wind-10 4.0 : 61.95% 53.32% 41.84% 30.30% 13.88% 1.08% 
Wind-10 8.0 : 41.68% 30.88% 19.01% 9.97% 2.09% 0.01% 
Wind-10 12.0 ؛‎ 28.04% 17.88% 8.64% 3.28% 0.32% 0.00% 
Wind-10 18.0 : 15.48% 7.88% 2.65% 0.62% 0.02% 0.00% 
Wind=10 21.0 : 11.50% 5.23% 1.46% 0.27% 0.00% 0.00% 
Wind-15 0.2 : 90.24% 89.55% 88.46% 87.01% 83.63% 73.42% 
Wind=-15 1.0: 83.25% 80.13% 75.33% 69.39% 56.89% 29.68% 
Wind-15 2.0 : 75.26% 69.73% 61.63% 52.29% 35.15% 9.57% 
Wind-15 4.0 : 61.51% 52.80% 41.25% 29.70% 13.42% 0.99% 
Wind-15 8.0 : 41.09% 30.28% 18.48% 9.58% 1.95% 0.01% 
Wind-15 12.0 ; 27.45% 17.36% 8.28X 3.09% 0.28% 0.00% 
Wind-15 18.0 : 14.99% 7.54% 2.48% 0.57% 0.028 0.00% 
Wind=15 21.0 : 11.08% 4.97% 1.36% 0.24% 0.00% 0.00% 
Wind-20 0.2 : 90.24% 89.55% 88.45% 87.00% 83.60% 73.37% 
Wind=20 1.0: 83.22% 80.09% 75.28% 69.33% 56.80% 29.57% 
Wind=20 2.0 , 75.21% 69.67% 61.55% 52.20% 35.04% 9.50% 
Wind=20 4.0 : 61.43% 52.71% 41.15% 29.59% 13.33% 0.98% 
Wind-20 8.0 : 40.99% 30.17% 18.39% 9.51% 1.93% 0.01% 
Wind=20 12.0 : 27.34% 17.27% 8.22% 3.06% 0.28% 0.00% 
Wind=20 18.0 : 14.90% 7.48% 2.45% 0.56% 0.02% 0.00% 
Wind=20 21.0 : 11.00% 4.92% 1.34% 0.24% 0.00% 0.00% 
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Tropical Frofile 
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TABLE OF A-COEFFICIENTS FOR SIX LOWTRAN PROFILES AT 3.5 um. 


Subarctic Summer Profile 


Visibility 40 AQ 20 15 10 

Wind = Û 0.900942 0.900967 0.900928 0.900916 9.900932 
Wind = & 0.900908 0.900991 9.900808 0.901004 0.897534 
Wind = 19 9.900920 0.900967 0.901247 4.900275 0.910749 
σπα = 15 0.900901. 0.900910 0.901037 0.900195 0.903370 
Wind = 20 0.901048 0.900953 0.900717 0.900709 0.889788 


Subarctic Winter Frofile 


Visibility 40 AO 20 15 10 

Wind = O 0296 5ذ‎ μωωωσσσο ٦۸960917 0.963549 30.955600 
Wind = 5 0.963548 0.965644 0.963652 9.963692 0.965209 
Wind = 190 0.96353594 0.963567 0.963535 0.963565 0.965684 
Wind = 15 0.96335 0.963605 0.963507 0.963501 0.956474 
Wind = 20 0.963387 0.963634 0.963474 9.963308 0.956391 

U.S. Standard Atmosphere Frofile 

Visibility 40 AQ 20 15 10 

Wind = 0 0.920811 0.920871 9.920850 0.920789 0.920864 
Wind = & 0.920841 0.920852 0.920746 0.920921 0.922365 
Wind = iù 0.920800 0.920765 0.920945 0.920259 0.919745 
Wind = 15 0.920802 0.920818 0.920675 0.920402 0.923206 
Wind = 20 0.920836 0.920866 0.920924 0.920379 0.914907 
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TABLE OF B-COEFFICIENTS FOR SIX LOWTRAN PROFILES AT 3.5 


um. 
Tropical Frofile 

Visibility 40 AO 20 15 10 
Wind = ο -0.078300 -0.079883 -0.083263 -0.086652 -0.093444 
Wind = 5 -Ο. 152564 —0.154864 —-0.199849 -0. 245554 -0.337061 
Wind = 10 -0.143678 -ū. 170376 -0.223794 -0.277764 -0.384270 
Wind = 15 -0.146815 -0.174519 -0.230264 -0.28697534 -0.401796 
Wind = 20 -0.148416 -0.176768 -0.233486 -0.291292 -0.404436 


Midlatitude Summer Profile 


Visibility 40 50 20 qu 10 

Wind - O —0. 0642753 —-0.065945 -0. 069515 -0. 072711 - ٦ 
Wind = 5 —0.118579 -0.14089? —-0. 185840 -0.251112 -Ο5Ε ٦ 
Wind - 10 -0.129858 -0.1563595 —0.209959 -0.263673 0.371240) 
Wind = 15 -0.132919 —0.160665 -0.216479 -0-Z273106.5 078 
Wind = 20 -0. 1.534394 -0.162714 -0.219772 -0.277832 -0.386820 


Midlatitude Winter Frofile 


Visibility 40 30 20 15 10 

Wind = 0 —-0.028635 -0.030311 -0.033691 -0.037106 -0.043922 
Wind = 5 -0.082925 -0.105237 -0.150164 -0.195458 -0.286223 
Wind = 10 -0.094293 -0.120931 -0.174555 -0.228508 -0.338264 
Wand = 15 -0.097416 -0.125235 -0.181256 -0.238064 -0.352103 
Wind = 20 -0.098874 -0.127235 —0.184428 -0.242007 -0.359376 
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TABLE OF B-COEFFICIENTS FOR SIX LOWTRAN PROFILES AT 


جھہٗ سے س سے سے سے س لے سے D D D rm mc me‏ س ل س mt sc ccm‏ کے لت ست اک سے کے کے سے سے سے س س س سے -e a ama cma ama Aee —XÀ‏ 


Wind 
Wind 
Wind 
Wind 
Wind 


Visibility 


— tr tere te te te tre ت ل‎ te te te tre tre cc ce crm tre توچ ین س سے سے سے سے سے سے لے ل للت للت کے اکا لے کے لے لے لے کے سے سل س س سل ل‎ eee eee ee ee ও ওল, ও এ 


سس سے سس س س س س کت سے لسا ل — — — — سے سے لے ل سے لست کے سے کس س لے سے لے لاسا ای me‏ اا ل ا سے سے سے لے لا لا لا للت لا لا کے لے کس سس س س রর‏ - ل س ل لے পার‏ ت لے کک 


—-0. 0488635 
-0 αυ 
-Ο. 114554 
πο 86 
ου 118852 


-0.017768 
-0. 071933 
-0. 083740 
-ο. 087009 
-0. 088548 


6 -0.0- 
1727ء بج 
z0.: 100272‏ 
102441 .0— 
103449 .0- 


— TABLE A.8 - 


Subarctic 


—-0. 050536 
—0.125554 
—0.140914 
—0. 1450953 
-0. 147091 


Subarctic 


-0.019474 
-0.094238 
70: 110310 
۲0 er) 50 
IIIS 


Summer Profi le 


—O. OS S881 
—0.170457 
—0.194472 
—Q. 2010072 
SBFF 


τ. 


-0.057260 
ο 77761 
-ο. 348144 
ডি) 
9 


Winter Frofile 


EOL 
—0.1.59144 
—0.164416 
mL 17/1452 
—0. 174743 


-0 . 026.575 
-0O. 184483 
—O. 2188698 
—-0. 228422 
—0. 2353146 


~ S 


— o u UM. 


-0.064021 
—0. 305416 
runc G 
2:025 50607 
—O. S 7220504 


SUDO 
-Ο. 274855 
TO SBT 
το. 340038 
το. 46 4 


U.S. Standard Atmosphere Frofile 


-0. 038290 
-0.004360 
-0.125884 
-0.128875 
-0.130276 
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=0: 041297 
—0.159 548 
ভিটা 7767 
πια 09 
—0. 184322 


hJ 


—O.04435 35 
03580 σοι 
—0.2294.59 
ο 9 9٦ 
το. ο 8687 


—-0. 050403 
—OQ. 296848 
πο 7 
—O0. 544040 


—Q. 545855 


APPENDIX B 


COMPUTER CODING TO INSTITUTE 
ATMOSPHERIC TRANSMITTANCE CALCULATION 


The following FORTRAN code is designed to institute the 
calculation of atmospheric transmittance on an incremental 
range movement basis. As the missile closes its target (s) 
periodic recalculation of acquisition probabilities based on 
updated range and atmospheric transmittance are necessary to 
realistically model an engagement. Listed here are the 
subroutines necessary to calculate the transmittance based 
on the value of the current range value in the model. 
Although these subroutines are written in FORTRAN, they are 
simple and easily convertable to other languages. Comments 


regarding user alternatives are contained therein. 


First thing to do is allow the user to select which 
Atmospheric Profile is desired. Then pass this to the 
ATMOSF subroutine to calculate specific coefficients. 
This input may be converted to disk-read input or 
NAMELIST input if interactiveness is not desired. 


Πρ endi im 


INTEGER ATFROF 


READ(111,*x) ATFROF 
111 FORMAT(’ Which Atmospheric Frofile is desired ?'// 
i is Tropical / 
Ze Midlatitude Summer 7 
ই Midlatitude Winter / 
4. Subarctic Summer ۶7 


30 Subarctic Winter / 
ó. U.S. Standard Atmosphere’ ,I1) 


ς- ς} 5 Οἱ Ι9 »- 


CALL ATMOSF CATFROF ) 


100 


- ——— پوو وى‎ > e 


οπου 093 


OOO QO 


1 


hJ 


The values of Wind and Visibility must already have 
been input to the program prior to the above CALL. 
ATMOSP calculates the ‘A’ and 'B' coefficients which 
are necessary to refer to the TRANSM function. 


== — —  — — — লহ س — — ”سے‎ লক ৮ 


SUBROUTINE ATMOSP (IT) 


DIMENSION A1 (6) ,A2Z(S) ,AS(6) ,A4(6) ,C1 (6) ,C2(5) , 
1 D1 (6) ,D2(5) , SIGMA (6) 


COMMON ACOEFF,BCOEFF ,VISIB,VWIND 
OPEN (UNIT=9,FILE='ATMPROF. DAT’ ,STATUS=‘OLD’) 


DO 1 I=1,6 
22210 05১৮) 21 01) 
READ(9,*) 7207) 
READ(9,*) AS(I) 
READ(9,*) A4(I) 
READ(9,*) 01 07) 
READ(9,*) 0201) 
READ(9,*) 01 01) 
READ(9,*) D2(I) 
READ(9,*) SIGMA (1) 

CONT INUE 


CLOSE (UNI T=9 ,STATUS=‘KEEP ' ) 


ACOEFF - SIGMA(CII) 

VCOEFF = A1(II) * A2(II) *«VIGIB * ASCIID *VISIB*x2Z. 
Hot ASCIID *VISGIB**3. 

WCOEFF = 1. + DICII)*#(1.- EXP(D2(1I1) #VWIND) ) 

IF (VWIND.LE.O.1) THEN 


VWCOEF=1.0 '/* No correlation if 
60 TO 2 '/* wind near zero. 
END IF 
VWCOEF = Ci(II) * VISIB ** C2(II) 
BCOEFF = VCOEFF * WCOEFF * VWCOEF 
RETURN 
END 


The two necessary coefficients (à or 60, B or @) are 
computed and global for referencing by TRANSM for 
transmittance calculation based on current range 
values. Now the function to perform equation 4.3. 
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FUNCTION TRANSM (RNG,AT) 


© The current range value is passed in as RNG, and AT is 
C the Atmospheric transmittance value passed back out. 


COMMON  RO,ACOEFF ,BCOEFF 
AT = ACOEFF * EXP (BCOEFF*RNG) 


RETURN 
END 


The file ‘ATMPROF.DAT’ referred to in subroutine ATMOSF 
above is an ASCII data file which contains all of the curve 
fitting coefficients listed in table 4.8. It 15 easily 
constructed with an editor, keeping in mind that if a name 
change is necessary, it must be performed in both places. 
The format is shown here and may be completed by referring 
to table 4.8. Ensure that no comments, blank lines or 


leading spaces occur when constructing the file. 


-0. 1391354 '/* Tropical Al 
6.2887E-O3 ! / ۴ a2 
—2.J44E-—O04 κ AS 
2-609 1E-06 te 74 
15. 24105 ! / +¥ eT 
-O. 716806 এ ER 
0. 895 κ 01 
۲٦ tx D2 
0.9383 κ 6 
-O. 1249832 '/* Midlat Summer Al 


etc. for others in the order specified. 


τος 
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